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Make speedier checks of recorders, control- 
lers and base or noble metal thermocouples 
in industrial] plants with the new three-dial 
8686 Portable Millivolt Potentiometer. Fea- 
tures such as a central reading window ... 
where measured values appear as a row of 
digits with a scale interpolation . .. simplify 
calibration of thermocouples and test meas- 
urements. The 8686 Potentiometer has: a 
wide operating range of —10.0 to +100.1 
mv and +1010 to +1020 mv for standard 
cell calibration; and a high accuracy of 
+(0.05% of reading +3uv) without refer- 
ence junction compensation, +(0.05% of 
reading + 6 nv) with ref. jet. comp. Write 


8686 Millivolt 
Potentiometer 


for Data Sheet E-33(1A). 


8692 and 8693 
Temperature 
Potentiometers 


If you want to make a variety 
of temperature measurments 
quickly with one flexible 
instrument, investigate the 
new time-saving 8692 Single- 
Range or 8693 Double-Range 
Temperature Potentiometers. 
Available in any of 24 inter- 
changeable temperature and 
millivolt ranges, these instru- 
ments read directly in degrees 
F or C on a scale 27%” long. 
Convenience features include: 
simplified range changes . . . 
only a screwdriver is needed 
to change a circuit panel, 
scale and binding post studs; 
automatic reference junction 
compensation . . . reference 
coil, built into circuit panel, 
compensates for thermocouple 
being used; accuracy ... 
+0.2% of range. Write for 
Data Sheet ND42-33(1A). 


instruments from the 


NEW 


CONVENIENCE-STYLED 
L&N LINE... 


Need a _ fast-operating, high-sensitivity, 
high-quality null indicator for use in re- 
search, testing and production checking? 
Here’s a new 9834 Guarded D-C Null Detec- 
tor having a short period of less than two 
seconds for source resistances up to 1000 
ohms, increasing to 4 seconds at 100,000 
ohms .. . ideal for measurements with 
guarded or unguarded potentiometers and 
bridges. Of rugged construction, this port- 
able, line-operated detector provides nu- 
merous convenience features which include 
four degrees of sensitivity, with a basic 
sensitivity of 0.2 nv/mm (0.3 pv/scale div.), 
and a noise level of less than +0.1 pv. 
Write for Data Sheet ED7(2). 


4992 Stenton Ave, 


LEEDS NORTHRUP 


Philedelphia 44, Pa 


8690 Millivolt 
Potentiometer 


Fast, accurate in-plant check- 
ing of thermocouples, record- 
ers and controllers is expe- 
dited by the new 8690 Port- 
able Millivolt Potentiometer. 
Design features such as “In- 
line Readout”, which permit 
measured values to be read in 
digits with a scale interpola- 
tion, accelerate temperature 
measurements. The 8690 pro- 
vides a wide range of —11.0 to 
+101.0 mv which eliminates 
input reversals; an accuracy 
of +(0.05% of reading + 20 
uv) without reference junc- 
tion compensation, + (0.05% 
of reading + 40 pv) with 
ref. jct. comp. Write for Data 
Sheet E-33(1B). 


9834 Guarded 
D-C Null Detector 
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The Franklin Medal (1914—Gold Medal) —This medal is awarded annually to those 
workers in physical science or technology, without regard to country, whose efforts, in the 
opinion of the Institute, acting through its Committee on Science and the Arts, have done 
most to advance a knowledge of physical science or its applications. 

The Elliott Cresson Medal (1848—Gold Medal).—This medal is awarded for discovery 
or original research, adding to the sum of human knowledge, irrespective of commercial value ; 
leading and practical utilizations of discovery; and invention, methods or products embody- 
ing substantial elements of leadership in their respective classes, or unusual skill or perfection 
in workmanship. 

The Howard N. Potts Medal (1906—Gold Medal).—This medal is awarded for distin- 
guished work in science or the arts; important development of previous basic discoveries, 
inventions or products of superior excellence or utilizing important principles. 

The John Price Wetherill Medal (1925—Silver Medal).—This medal is awarded for 
discovery or invention in the physical sciences or for new and important combinations of 
principles or methods already known. 

The Edward Longstreth Medal (1890—Silver Medal).—This medal is awarded for 
inventions of high order and for particularly meritorious improvements and developments in 
machines and mechanical processes. In the event of an accumulation of the fund for medals 
beyond the sum of one hundred dollars, it is competent for the Committee on Science and the 
Arts to offer from such surplus a money premium for some special work on any mechanical 
or scientific subject that is considered of sufficient importance. 

The Louis E. Levy Medal (1923—Gold Medal).—This medal is awarded to the author 
ofa paper of especial merit, published in the JouRNAL OF Tue Franxun Institute, prefer- 
ence being given to one describing the author’s experimental and theoretical researches in a 
subject of fundamental importance. 

The George R. Henderson Medal (1924—Gold Medal).—This medal is awarded for 
meritorious inventions or discoveries in the field of Railway Engineering. 

The Walton Clark Medal (1926—Gold Medal).—This medal is awarded to the “author 
of the most notable advance in knowledge or improvement in apparatus, or in method con- 
cerning the science or the art of gas manufacture or distribution or utilization in the produc- 
tion of illumination, or of heat, or of power.” 

The Frank P. Brown Medal (1938—Silver Medal).—This medal is awarded to in- 
ventors for discoveries and inventions involving meritorious improvements in the building and 
allied industries. 

The Newcomen Medal (1943—Gold Medal).—This medal is awarded, not oftener than 
once in three years, for achievement in the field of Steam. 

The Francis J. Clamer Medal (1943—Silver Medal).—This medal is awarded at least 
ence in five years for meritorious achievement in the field of Metallurgy. 

The Stuart Ballantine Medal (1946—Gold Medal).—This medal is to be awarded in 
recognition of outstanding achievement in the fields of Communication and Reconnaissance 
which employ electromagnetic radiation. 

The Boyden Premium (1859).—This premium is awarded not oftener than once in five 
years to any resident of North America who has recently made a notable experimental deter- 
mination of the speed, in free space, of radiation in any region of the entire spectrum. 


The William M. Vermilye Medal (1937—Bronze Medal).—This medal is awarded not 
med than biennially in recognition of outstanding contribution in the field of Industrial 
anagement. 


_ . The Certificate of Merit (1882).—A Certificate of Merit is awarded to persons ad- 
judged worthy thereof for meritorious inventions, discoveries or improvements in physical 
processes or devices. 

For further information relating to these awards apply to the Secretary, Committee on Science and the Arts 
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Announcing publication of 


The Proceedings of the Symposium on 
—GAS-COOLED REACTORS-— 


held February 10-11, 1960 


sponsored jointly by The Franklin Institute and the Delaware 
Valley Section of the American Nuclear Society. 


Twenty-one papers are featured by well known 
U. S. research people and several foreign scientists 
who have reviewed the developments in this new type 
of atomic power plant. The symposium aimed to 
clarify the problems which must be solved before gas 
cooling can be used to its fullest potential. 


Sir William Cook, C.B., of the U. K. Atomic Energy 
Authority reviewed the British program in gas-cooled 
reactors; a representative from the Swiss firm, Brown, 
Boveri & Cie., discussed their work on a pebble bed 
reactor which they are —_— for the West German 
Government. Representatives of the Australian Atomic 
Energy Commission and from West Germany pre- 
sented programs being carried out in their respective 
countries. Discussions are included at the end of 


most papers. 


The sessions were devoted to the following topics: 
| Status of Major Gas-Cooled Reactor Programs 
Il Fuel Elements and Materials 
111 Components and Loop Studies 


Issued as Monograph No. 7 under the auspices of the JOURNAL OF THE FRANK- 
LIN INSTITUTE, the Proceedings are now available. Send your order (at $5.00 per 


copy) to 


JOURNAL OF THE FRANKLIN INSTITUTE 


20th & Parkway 
Philadelphia 3, Pa. 
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there is 
no good 

and original 
observation” 
— Charles Darwin, naturalist 


Man’s search for scientific knowledge and understanding. has its 
| taproots in the above thought expressed by Darwin in a letter 
-_ tohis distinguished contemporary, Alfred Russel Wallace, in 1857, 


Speculation —intuitive contemplation guided by past discoy- 
' > eries—led Darwin to his famous observations set forth in Origin 
of Species. Similarly, it led Alexander Graham Bell to the inyen- 
| tion of the telephone —and has since led to dozens of major ad- 
vances in the field of electrical communications. 


a At Bell Telephone Laboratories, the puzzling flow of current 
~~ im semiconductors provoked speculation which yielded the tran- 
sistor—and a Nobe! Prize. Speculation about the behavior of the 
electron led to experimental proof of its wave nature — and another 
_ Nobel Prize. “Brains” capable of guiding missiles and spate 
pfohesicame first from speculations at Bell Laboratories. 


2 Today. Bell Laboratories scientists and engineers are more 
_. keenly aware than ever of the importance of speculative thinking. 
The far-reaching scientific and technological developments of to- 
morrow are already the subject of advanced research, Amongthem 
are radically new materials and devices—basically new switching 
systems, transmission via satellites, and waveguide networks able 
to carry hundreds of thousands of voices simultaneously. 


Through informed speculation about Nature’s laws, Bell Lab- 
oratories will continue to search for the “good and original obser- 
vations” so vital to the ever-improving Bell Telephone System, 
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NEW GYROS FOR OUR SUBMARINES * 


BY 


FREDERICK D. BRADDON ! 


ABSTRACT 


This report describes the present day use of gyro compasses on naval vessels. 
The nuclear powered submarine with its capability for long sustained submergence 
and long range cruising, presents new and challenging problems of navigation. Two 
new types of gyro compass have been introduced by the Navy in recent years to meet 
the increasing navigational and ordnance requirements of the fleet. Several hundred 
installations of these equipments have been made over the past four years in both 
submarines and surface ships. This background of experience serves to establish 
thoroughly the operating characteristics of these equipments. This talk describes 
the design features, operating principles, and performance capabilities of the Mark 19 
Gyro Compass and its companion instrument, the Mark 23 Gyro Compass. Also 
described are some of the specialized problems of polar navigation and the techniques 
employed by the U.S.S. Nautilus, the Skate, and the Sargo in the use of their gyro 
compasses for their successful polar transits. The Mark 19 Gyro Compass represents 
the first fully stabilized directional reference with automatic north seeking control, 
and is the direct forerunner of the Ship Inertial Navigation Systems. The Mark 23 
Gyro Compass is unstabilized and indicates only azimuth. It serves for smaller 
vessels and as an auxiliary compass in combination with the Mark 19. I will try to 
show the continuing function of the gyro compass in relation to this more complex 
type of GYRO NAVIGATOR™ now being employed on the new Polaris submarines. 


On every naval vessel one will find an instrument known affection- 
ately as the Gyro. This is the compass on which the ship depends to 
guide it on its way. Each man aboard knows that his safety and his 
mission depend upon this Gyro every minute of the day. It functions 
unattended and incessantly for months on end. It is the central datum 
to which the ship is steered, and to which all navigation, surveillance 
and ordnance functions are referred. 


* Lecture delivered May 4, 1960, at a Joint Meeting of the Philadelphia Science Teachers 
Association and The Franklin Institute. 
! Chief Engineer, Marine Division, Sperry Gyroscope Company, Syosset, N. Y. 


(Note—The Franklin Institute is not responsible for the statements and opinions advanced by contributors in 
the Journat.) 
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The Gyro is not a conspicuous instrument, for it is small and quiet 
and is generally installed in a protected area deep within the vessel. 
It requires only electrical connections for its operation. Its output 
information is transmitted electrically to the remote navigation, sur- 
veillance and ordnance stations which need to know true-north and the 
zenith. This equipment is technically designated as the ship’s Master 
Gyro Compass. 

Many forms of gyro compass have been developed and employed 
by the Navy since the first installation on the U.S.S. Delaware in 1911. 
Changing requirements and changing technology have resulted in a 
succession of designs that have progressively led to improved accuracies, 
size reduction and a high degree of operating reliability. It is the 
purpose of this talk to describe the two currently manufactured Navy 
type compasses known as the Mark 19 and Mark 23 which are in wide 
operational use. Figure 1 shows the design and relative size of these 
equipments. 


Fic. 1. Composite photos of MK 19 and MK 23 gyro compasses. 


The origin of the Mark 19 dates back to World War II when we were 
confronted with the problem of precisely stabilizing and directing a high 
resolution shipborne radar designed for aircraft surveillance and air 
traffic control on carriers. We were amazed to discover the inadequacy 
of existing shipboard gyro compasses and stable elements for this 
purpose. Gyro compass accuracy was reasonably adequate for navi- 
gation, but transient perturbations in the output data hindered its use 
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for azimuth stabilization. The stable elements lacked vertical ac- 
curacy due to instrumental errors caused by gyro friction, mass shift 
and inadequate compensation of maneuvering and seaway accelera- 
tions. A comprehensive study was initiated to analyze this problem 
from which evolved the concept of a unitary three-axis gyroscopic 
element to replace the two separate instruments for azimuth and vertical 
indication. 

Information was collected on the various forms of gyro compass 
and gyro vertical which were in general use in 1944, or which had been 
previously developed to the point of successful operation. We cata- 
logued the detailed design and operating characteristics of these 
equipments under the headings listed in Fig. 2. From these data we 


PHYSICAL PROPERTIES 
NUMBER OF GYROS STABLE 
ANGULAR MOMENTUM 
EQUIPMENT SIZE 
EQUIPMENT WEIGHT Stan! 


DESIGN FEATURES 


GYRO SUSPENSION 
GRAVITY CONTROL C 
DAMPING OF OSCILLATIONS 0 
PLATFORM STABILIZATION Y am 

DATA TRANSMISSION R a7 

LATITUDE COMPENSATION 0 Ve 


PERFORMANCE CHARACTERISTICS 
INTERCARDINAL ROLLING ERROR 
SPEED ERROR COMPENSATION 0 ow 
TURN ERROR COMPENSATION 


GIMBAL ERRORS 

DATA TRANSMISSION ERRORS P 

DAMPING-LATITUDE ERROR EFFECTS A. wae 
OPERATING FEATURES $ 

EASE OF MAINTENANCE E 

OPERATING DEPENDABILITY 

SERVICE LIFE 


MANUAL SUPERVISION REQUIREMENTS 


ore 


UNDER STATIC CONDITIONS 
UNDER DYNAMIC CONDITIONS 
IN POLAR REGIONS 
FOR NAVIGATION REQUIREMENTS 
FOR ORDNANCE REQUIREMENTS 
- — 


Fic. 2. Evaluation of gyro compasses and stable vertical equipment as of 
November, 1944. 


attempted to analyze and evaluate the theoretical and practical aspects 
of each design feature. On this basis we established the initial specifica- 
tions for what was to become the MK 19 Gyro Compass. This tabula- 
tion further indicated that all the gyro compasses used gyros with a 
horizontal spin axis. Figure 3 shows that Sperry and Brown used a 
single gyro in the meridian plane, while Anschutz and Arma used two 
horizontal gyros linked to subtend the meridian plane. All the gyro 
verticals (sometimes referred to as stable elements) used a single gyro 
with a vertical spin axis. At this point a very significant observation 
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82 FrepericK D. Brappon J. 
was made. The most precise level reference on a ship is the level attitude 
of the spin axis of the gyro compass. Because of the extreme freedom 
of suspension of the gyro compass about the vertical axis, the gyro is 
free of perturbations in level. The gyro compass was therefore capable 
of more accurate indication of tilt in the meridian plane than any of the 
gyro verticals. From this observation it was deduced that an identical 
gyro with horizontal spin axis slaved 90° from the meridian gyro, as 
shown in Fig. 4, would sense tilt in the E-W plane with equal precision. 


ANSCHUTZ AND ARMA SPERRY AND BROWN GYROCOMPASSES 
GYROCOMPASSES (WITH HORIZONTAL SPIN AXIS) 


ALL GYRO VERTICALS 
(VERTICAL SPIN AXIS) 


OF EARTH 


EQUATOR 


Fic. 3. Gyro arrangements for azimuth and vertical determination. 


Two intersecting horizontal lines thus define the horizontal plane to 
establish a zenith reference and the north seeking control of the meridian 
gyro continued to point its spin axis to true north. The MK 19 con- 
cept thus evolved as a precise and compact single instrument to replace 
the separate indicators. We originally called this instrument the 
“Zenith Meridian Indicator,’ and later conformed to Navy terminology 
using the designation ‘““Gyro Compass Equipment Stabilized MK 19.” 
This title has tended to de-emphasize the vertical indicating capability 
of the instrument. 

As the result of this study and extensive experimentation and proto- 
type manufacture, the present MK 19 Mod 3 design was established 
having the basic features listed in Fig. 5. The idea of the composite 
vertical and azimuth indicator, the use of electrical methods of applying 
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hic. 4. Concept of unitary three-axis stable element. 
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gravity coupling, and ship speed compensation by gyro torquing were 
all radical innovations. 

Before describing the design characteristics of this compass, it is 
well to review the elementary principle of operation of the gyro com- 
pass that slaves it to the Earth and gives it the ability to automatically 
align itself to true north. A gyro compass serves to indicate the local 
meridian. Referring to Fig. 6, the meridian plane is defined by two 
intersecting lines: (1) the axis of rotation of the Earth, and (2) the 
plumb-line at the ship’s location. We use a gyroscope to sense the 
rotation of the Earth. On this gyro we mount a pendulum, liquid in 
containers, an accelerometer, or a spirit level to sense the attitude of the 
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Fic. 6. 


North-seeking principle of the gyro compass. 


gyro to the gravity vector. As the Earth turns under the gyro, which 
tends to remain fixed in reference to the stars, the spirit level detects a 
tilt. Inthe MK 19 and MK 23, an electrical signal is derived from this 
level and an electrical torque applied to direct the gyro toward the 
meridian. The gyro thus nulls on the horizontal component of Earth 
rotation to settle in the meridian plane. A natural error in azimuth 
indication results from the northerly velocity of the ship. The instru- 
ment can sense only the total angular velocity of the Earth and ship. 
Figure 7 shows this effect to cause the gyro to settle to the west of the 
meridian for northerly transits. In all previous types of gyro compass, 
the gyro spin axis was deflected from the meridian by a relatively large 
amount as a function of the ship’s speed, heading, and latitude. This 
error was computed and corrected in the output data. As a result, 
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perturbations in velocity log data caused instantaneous transients in the 
output data. In the MK 19 and MK 23 equipments, a vertical torque 
proportional to northerly speed is applied electrically to tilt the gyro 
at a rate equal to ship’s northerly velocity. This control completely 
compensates the northerly speed error, smooths the effect of rough 
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Fic. 7. Northerly speed azimuth error. 


| 


LEVELING 
\. CONTROL 


< 


~*~ AZIMUTH 
CONTROL 
AMPL 


SLAVE GYRO TILT SIGNAL 


TYAS 


SLAVE GYRO 
LEVELING 
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velocity data, and maintains the gyroscopic element fixed in the 
meridian plane. 

Without too much detail, Fig. 8 shows schematically the composition 
of the MK 19 Stable Element. The level signal of the meridian gyro is 
amplified and applies a horizontal torque to direct the system to the 
meridian as previously described. Simultaneously, the same signal 
applies a smaller torque about the vertical axis to level the spin axis. 
This represents a damping function in the servo loop. The Slave Gyro 
level applies a similar leveling torque to its gyro while the vertical axis 
pickoff applies a horizontal torque to maintain the alignment of this 
gyro normal to the meridian gyro. 


AZIMUTH, 
FOLLOWUP / 
AMPL 


\ 


- 


MERIDIAN 
GYRO 


AZIMUTH 
FOLLOWUP 
PICK-OFF 


SYNCHRO TRANSMITTER 
(HEADING DATA) 


AZIMUTH PHANTOM 


TAY 
\\ 


SL AVE \ 
GYRO 
Vy 


Fic. 9. Azimuth phantom servo drive 


Figure 9 shows the addition of the azimuth servo which maintains 
the phantom in alignment with the meridian gyro. 

Figure 10 shows the addition of the roll and pitch gimbals with their 
corresponding servo drives and deck tilt data transmitters. The con- 
trol signals for roll and pitch stabilization are derived from the hori- 
zontal pickoffs on the two gyros. These signals are necessarily trans- 
formed from terms of E-W tilt and N-S tilt to the corresponding roll 
and pitch values by use of an electrical resolver on the azimuth drive. 

Figures 11-15 show the construction features of the MK 19 Equip- 
ment. The rotor shown in Fig. 11 is about 5} inches in diameter, 
weighs 8 Ib., spins at about 400 revolutions per second, and represents 
an angular momentum of about 90 X 10° gm cm?/sec. The design 
stresses absolute symmetry of construction, high electrical efficiency, 
and an efficient lubrication system. 
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;. 10. Platform gimbal and stabilization servos. 


Fic. 11. Gyro motor construction. 


The rotor is sealed within a spherical case, shown in Fig. 12. This 
assembly is accurately balanced under fluid for neutral buoyancy and 
to bring the centroid of flotation coincident with the center of mass. 
Many of these units have operated without benefit of maintenance for 
periods exceeding 12,000 hours of service. An interesting comparison 
is to consider 12,000 hours of driving your car at 30 mph., thus travelling 
360,000 miles, without an oil change or adjustment of any kind! This 
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Fic. 12. Gyro float assembly. 
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Fic. 13. Gyro gimbal and level assembly. 


12,000 hours (or almost 1} years) of continuous operation at 400 
revolutions per second represents 17 billion revolutions of the rotor with 
a peripheral travel of the rim of 3} million miles. One of the early 
installations on an ice breaker logged over 19,000 hours of satisfactory 
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| 


operation. Two design features contributing to this reliability of 
operation provide a thermally controlled lubrication system and positive 
control of the rotor bearing preload. 

The sealed gyro unit is then installed in a sensitive gimbal shown in 
Fig. 13. The electrolytic level, angular pickoff and slip rings are then 
added to this member. This assembly then mounts on the base of the 
oil filled enclosure as pictured in Fig. 14. The Slave Gyro is similar in 
all construction details. 

These completed units are fully adjusted and tested separate from 
the system, and are stored and installed as interchangeable assemblies. 
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14. Meridian gyro mounting. 


Figure 15 shows the meridian and slave units installed on the 
stabilized platform ready for system tests. (Note that the shock 
mounts are incorporated within the binnacle and a mechanical panto- 
graph linkage serves to provide positive alignment of the stabiliza- 
tion axes.) 

The very slow response characteristics inherent in the gyro control 
system make it mandatory to avoid any interruption in operation of the 
equipment. Even a momentary interruption will destroy system 
accuracies for hours to follow. Provision is made in equipment design 
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to guard against this event. Figure 16 shows the construction of the 
Control Cabinet of the MK 19, which houses the servo amplifiers, 
power converters, computers, electro-mechanical integrators, indicators, 
and test circuits. Magnetic amplifiers are used for reliability in the 
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Fic. 17. MK 23 gyro compass gyro and gimbal arrangement. 
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Fic. 19. MK 23 master gyro compass 


critical stabilizer servos. Electronic amplifiers are applied in parallel 
to guard against component failure. Monitor and checkout circuits 
permit failure or replacement of any component without interruption 
to system operation. 

The MK 23 Gyro Compass, represented in Fig. 17, incorporates all 
of the features of the MK 19 meridian gyro, but substitutes a simple 
pendulous gimbal suspension for the roll and pitch servo stabilized 
platform. With other design compromises, the equipment represents 
about one fourth the complexity, size and cost of the MK 19 Gyro 
Compass. It, of course, provides only azimuth data. The schematic 
shows the same control features previously described for the MK 19. 

Figures 18, 19, and 20 show the construction details of the MK 23 
Gyro Compass. Its operating reliability is favored by the avoidance 
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of the critical stabilizing servos, by its over-all simplicity, and by its 
low demand on ship’s electrical power. When used as an emergency 
back-up for MK 19 installations, care should be exercised to see that 
the MK 23 operates from an independent source of ship's power. 
Interruptions in the switching of ship’s power is the most frequent cause 
of error in compass operations. 

The generally accepted performance accuracies of the MK 19 and 
MK 23 Gyro Compasses exceed the design specifications by an appre- 
ciable amount. Figure 21 represents a summary tabulation showing 
the static and dynamic performance capabilities of these equipments. 
To generalize, we might classify the MK 19 as accurate in azimuth to 
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Fic. 21. MK 19 and MK 23 gyro compass performance accuracy (at 40° latitude). 


one tenth degree when operating in mid latitudes under normal ship- 
board conditions. The corresponding accuracy of the MK 23 is about 
one half degree. The .05° azimuth repeatability of the MK 19 is a 
significant figure in establishing the quality of the instrument. This 
value indicates that the gyro is consistently, day after day, measuring 
Earth rates of 11.7°/hour X sin .05°, or .01°/hour. 

To the gyro engineer and navigator, the performance of the MK 19 
and 23 Gyro Compasses on the Nautilus polar transit was almost as 
spectacular as the mission itself! The actual compass log entries 
show that on August 1, 1958 at 5:00 p.m. and latitude 71—51.8, the sub- 
merged ship was pointed to the North Pole. With a MK 19 heading of 
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000 and in automatic steering, the ship cruised for seven hours that day, 
all the next day, and for 23} hours the following day, to latitude 88.46 
without alteration of course. With the Pole indicated as less than 90 
miles ahead, a small course correction of 3° was made to correct an 
indicated offset of only 43 miles after 1100 miles of transit! It is 
equally impressive to note that the divergence of the MK 23 from the 
MK 19 never exceeded +1°. This circumstance was especially reas- 
suring to the operation. 

The pictorial representation of this voyage, shown in Fig. 22, (and 
with acknowledgment to National Geographic) shows that after a total 
submerged transit of over 96 hours and 1,800 miles, the initial celestial 
position fix checked within 10 miles of the dead reckoning position. 


\ 


\\, 


Fic. 22. Nautilus polar transit. 


The Nautilus expedition was planned in a manner to favor gyro 
compass operation with straight courses, no maneuvering accelerations, 
and minimum time of transit. MK 19 performance on the Skate and 
Sargo was equally satisfactory below 80° latitude, but showed fairly 
high drift rates while maneuvering for days in the extreme polar region. 

Although the MK 19 Gyro Compass was not designed to be operated 
as a Directional Gyro for grid-type navigation in the polar region, the 
equipments on the Nautilus, Skate and Sargo have been modified for such 
use. The present design is not ideally suited to this type operation 
because of the mass shift effects of the bubble in the electrolytic level 
and the variable sensitivity of this device. A Navy sponsored develop- 
ment program now in progress will greatly increase the Directional 
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Gyro stability of the MK 19 for polar navigation. Figure 23 illustrates 
the bubble level presently used and the force feedback pendulum which 
will replace the spirit level. In the spirit level, the variable wetting of 
the electrode, with tilt, alters the balance of the resistance bridge circuit. 
A clean signal with high tilt sensitivity is provided by a small and simple, 
reliable device. Thermal sensitivity, shock sensitivity, mass shift 
effects, and response lag are, however, considered to limit system per- 
formance accuracies, particularly in the Directional Gyro mode of 
operation. 
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Fic. 23. Liquid level versus torque feedback accelerometer. 


The force feedback pendulum is somewhat more complex, but is 
expected to avoid all of the operating limitations of the spirit level and 
improves both the static and dynamic accuracies of the MK 19 System. 
This modification is equally applicable to the MK 23 System. Other 
features in this redesign of the MK 19 Gyro Compass will further sim- 
plify the system and lead to improved performance accuracies. 

The requirement for precise zenith and azimuth information, with 
absolute accuracy and stability for satisfying the increasing demands of 
ordnance and navigation, represents a continuing need within the Navy 
for the gyro compass. In this talk I want to help you to fully under- 
stand those features that make the gyro compass ideally suited to the 
normal navigational needs of the nuclear submarine. 
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The gyro compass is automatic in its alignment to the meridian 
plane through the Earth’s polar axis. It requires only to feel the force 
of gravity and to sense the daily rotation of the Earth about its polar 
axis. These two natural quantities persist in all parts of the Earth and 
for all time. They exist inside the vessel as well as outside, and pene- 
trate the depths of the ocean. They are unaffected by daylight or 
nightfall, or by the seasons of the year. They cannot be altered or 
obscured by weather or by man-made interference. The gyro compass 
thus provides our ships with secure, continuous and precise compass 
information independent of the stars, magnetic fields-or radio trans- 
mission. For compass information, a submarine need never break the 
ocean surface and thereby expose its presence to the enemy. 
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Navigational accuracy of dead reckoning is primarily limited by the 
lack of means to measure speed precisely over the bottom. A ship's 
heading error of 0.1° represents an offset from course of less than 2 miles 
in position after 1000 miles of transit. If we could measure ground- 
speed precisely, this 2 miles would be our position error. Today with 
uncharted ocean currents and limitations in measuring speed through 
the water, an uncertainty of 1 knot in estimating bottom speed can 
develop position errors up to 30 miles or more in a 1000-mile mission. 
Some day we will find a secure means to measure ground speed. 

The MK 19 made available for the first time in a submarine a pre- 
cise zenith reference. A logical and relatively simple outgrowth of 
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this program was to bring this information into the periscope so that 
star altitude measurements could be made for a celestial position fix. 
This equipment constitutes a periscopic sextant known as the Sub- 
marine Celestial Altitude Recorder (SCAR), shown in Fig. 24. This 
equipment is installed on a number of our advanced submarines and for 
the first time permits the vessel to establish its position by celestial 
reference without surfacing. 

To place the MK 19 and MK 23 program in perspective I have 
tabulated a few of the mileposts in the long history of gyro instrumenta- 
tion. We are told that the unique characteristics of the gyroscope were 
first analyzed by Newton, based on observations and theories of astron- 
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Fic. 25. Chronology of gyro navigation. 


omers relating to the cyclic precession of the equinoxes. Figure 25 
shows that Bohenberger in 1827 is credited with the conception of the 
gyroscope as a mechanical device. Lang, about 1846, suggested its 
possible use as a gyro compass. It remained for Foucault in 1852 to 
demonstrate this effect and expound the theory of the device. Trouve, 
Kelvin, van den Bos and Siemens contributed significant experimenta- 
tion from 1865 to 1885. Martienssen in 1905 developed, for the first 
time, an accurate mathematical analysis of the gyro compass. About 
five years later Anschutz in Germany and Elmer Sperry in New York 
produced the first practical shipboard equipments. 

Following the first Navy gyro compass installation in March, 1913 
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the Navy has continued to purchase various types of navigational and 
ordnance gyro compasses. From the early twenties until recent years, 
Arma was an alternate supplier to the Navy. Figure 25 shows that 
only MK 14, 19, and 22 types are in current production. The approxi- 
mate number of Navy ship installations for each type is shown. The 
total over a span of 50 years represents something over 10,000 
equipments. 

The very complex Ship Jnertial Navigation System (SINS), as 
now being produced for our new fleet of Polaris submarines, represents 
a new phase in ship navigation. This equipment incorporates the geo- 
graphic directional reference system of the stabilized gyro compass 
producing precise zenith and azimuth data. In addition, it provides a 
continuous record of ship position in terms of latitude and longitude, 
from which is derived a measure of true ground speed and other func- 
tions. Ship position change is detected by self-contained sensors which 
respond to the prevailing accelerations as measured in the meridian and 
East-West planes. Phenomenal precision is required in all of the out- 
puts of this instrument to satisfy the specialized ordnance use for which 
it serves in directing the Polaris intermediate range missile. In its 
present concept, the GYRO NAVIGATOR as used in SINS is neces- 
sarily monitored by an array of support equipment and requires con- 
tinuous supervision which can only be justified for the most vital 
missions. Most of our Navy will continue to depend on the gyro 
compass to meet its normal navigation, surveillance, and ordnance needs. 
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HARMONIC AND SUBHARMONIC RESPONSE OF AN ON-OFF 
CONTROL SYSTEM TO SINUSOIDAL INPUTS * 


BY 


B. A. FLEISHMAN ' 


ABSTRACT 


We derive here periodic solutions, of harmonic and subharmonic type, of a 
differential equation governing the operation of a relay servomechanism which is 
subjected to a sinusoidal input yo sin wt. Each solution is a “formal” solution in the 


sense that it is constructed under an a priori assumption about the manner in which 
it changes sign, and this assumption can be verified only after one has the solution. 
A significant part of the analysis, then, is devoted to the determination of suitable 
conditions on the parameters yo and w which insure the validity of the formal solution. 


INTRODUCTION 
Consider the second-order piecewise-linear ordinary differential 
equation 
¢+x+sgnx = — jj (1) 


where - = d/dt, -- = d*/dt*, and sgn x is defined by 


(x > 0) 


x 

sgn x = signum x = > = 
x | 
\ 


—] (x < 


The quantities x, y and ¢ are assumed to be dimensionless. 

It is shown in Appendix A that the differential equation of a simple 
servomechanism, with proportional plus discontinuous control, may 
be reduced to the form (1); ¢ represents the time variable, y = y(t) the 
(known) input function, and x = x(t) the (unknown) error function, 
that is, deviation of output from input. 

The basic analytical problem associated with any control system is 
to predict the error x(t) (or, equivalently, the output 2(4) = x(t) + y(t) 
as a function of time, when the input function y(t) and the initial error 
and error rate are prescribed. Because it is not linear, Eq. 1 cannot be 
solved exactly for arbitrary inputs y(t). Most analyses of nonlinear 
control systems deal with step (y = constant) and ramp (dy/dt = con- 
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stant) inputs, principally by phase-plane techniques (for example, see 
(1)*). Except by numerical means, very little has been done with other 
continuously varying inputs y(¢), particularly for differential equations 
with discontinuous nonlinearities, where such classical techniques as 
iteration and perturbation are inapproprate. 

In this paper sinusoidal inputs are considered. Even so, because 
(1) is not linear, we cannot find an exact solution for arbitrary initial 
data. Rather, we study only periodic solutions. It is shown that 
when y(t) is a sine function, the system governed by (1) is capable of 
steady-state periodic response, of harmonic or subharmonic type, de- 
pending upon the initial conditions and the values of the parameters 
in the input function. A harmonic solution is one with the same period 
as the input, while a subharmonic response is one whose smallest period 
is an integral multiple of that of the input. These periodic solutions are 
obtained by exploiting the piecewise-linear nature of Eq. 1. 

The occurrence of subharmonic oscillations is a significant feature 
of nonlinear systems not shared by linear ones. If its amplitude is 
large, such an oscillation may have destructive effects. Several 
writers have considered this problem. West and Douce (2) have 
analyzed subharmonic responses in nonlinear control systems by means 
of the describing-function method. Very recently Ogata (3) has 
studied the phenomenon analytically and experimentally (by use of an 
analog computer) in a group of systems in which the output of the non- 
linear element feeds into linear elements. 

The analysis presented here deals with differential equation (4) 
below; in this system the discontinuous and proportional control 
elements are in parallel, that is, they act simultaneously on the error 
signal x(t). We are able to derive explicitly, and verify, the periodic 
solutions, pointing the way to a study of periodic responses to general 
periodic (not necessarily sinusoidal) inputs. This extension will be 
carried out in a future investigation. One of the principal results of this 
paper is the exact analytical form, in the case of the control system 
studied here, of the subharmonic oscillations. 

Since we only consider inputs of the form 


(3) 


y(t) = yosin ot, 


where the amplitude yo and frequency w are constant, Eq. 1 may be 
written 


#+x+ sgnx = yo’ sin wf. (4) 


Clearly, once we find any particular solution x(/) of Eq. 4, representing 
the error as a function of time, we shall also have the output, z = x + y 
=x + yosin otf. 


? The boldface numbers in parentheses refer to the references appended to this paper. 
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The harmonic case is treated in Part I. First a ‘“‘formal’’ steady- 
state solution of Eq. 4 with period 27/w is constructed, then restrictions 
are imposed on the parameters yo and w to insure that this periodic 
function is indeed a solution of the differential equation. In Part II, 
by suitable modification of the harmonic solution, we obtain a formal 
subharmonic solution (of period 2mx/w, m an integer), and again we 
place restrictions on the parameters to insure the validity of this formal 
solution. The principal results are contained in theorems at the end of 
Parts I and II, and in Part III, “Concluding Remarks.” 


I. HARMONIC CASE 


1. Derivation of Formal Solution 


A harmonic solution of Eq. 4 can be constructed as follows. _ If two 
functions x,(t) and x,(t) satisfy the respective equations 


+ = Yow" sin wl (5) 


and 


+ sgn xy 


then upon adding these equations we see that the function 


x(t) = x,/(t) + x,(t) 


is a solution of Eq. 4, provided that sgn x, = sgn x,, that is, provided 
that x, and x, always have the same sign. 

Since there is no damping in the system, let us suppose that x,, x; 
and x, are all in phase (or 180° out of phase) with sin wf. The particular 
solution of Eq. 5 with this property is just 


x,;(t) = sin wt. (8) 


As for x, and x,, they must vanish at ft = kr/w, k = 0,1,2,---. Fur- 
thermore, since they are to have period 27/w and be in phase with sin of, 
it is not unreasonable to assume that they maintain one sign during the 
half-cycle 0 < t < r/w, the opposite sign during the next half-cycle, 
etc. To be specific, we assume x, (and therefore x,) > 0 during the 
first half-cycle. 

Thus, in view of Eq. 6 we seek a periodic solution x,(t) of the 
equations 


#+x=-1 (x > 0) (9a) 


#+x=+1 (x < 0) (9b) 
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having the period 2r/w and satisfying 
x(0) 0), z(0) = Vo > 0, (10) 


where the constant v» remains to be determined. We can gain our 
objective most elegantly by way of a geometrical discussion. 

Multiplying both sides of Eq. 9a by z and integrating with respect 
to t, we get 


z? + (x + 1)? = constant. (11a) 
In other words, all trajectories of Eq. 9a in an x, #-plane (phase plane) 
are arcs of circles centered at x = — 1, = 0;if x(0) = x, and = wo, 


the solution of (9a) with these initial values is represented by the circular 


x 


Fic. 1. Construction in phase plane of solution of Eqs. 9 with period 2x/w. 


arc centered at (—1, 0) which passes through the point (x, vo). When 
the positive x- and #-axes extend, respectively, upward and to the right 
(as in Fig. 1), an arc is traversed in the counterclockwise direction about 
its center as ¢ increases, for when # = dx/dt is positive x must increase 
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with increasing ¢. Similarly, all phase-plane trajectories of Eq. 9) are 
arcs of circles 
#2 + (x — 1)* = constant (115) 


centered at x = +1, = 0, again always traversed in the counter- 
clockwise direction. 

Now note that Eqs. 9 contain restrictions on the sign of x: Eq. 9a 
applies only when x is positive, Eq. 95 when it is negative. Thus a 
circular arc of type (11a) may be a portion of a trajectory only while the 
trajectory is in the upper half-plane. When such an arc, as part of a 
trajectory, intersects the z-axis the trajectory continues beyond this 
point of intersection as an arc of type (115), until this arc in turn inter- 
sects the z-axis, etc. 

A final remark on the nature of all trajectories of Eqs. 9 concerns 
the time it takes for a given circular arc to be traversed. If Eqs. 11 
are written in the form 


[const. — (x + 1)*}! 
and the variables are separated, we see immediately upon integrating 
that the time elapsed during the traversal of any (circular arc) trajec- 
tory is just equal to the angle subtended by the arc at its center. 

We are now in a position to construct the required solution x,(¢) of 
Eqs. 9. First of all, because of the initial conditions (10), the tra- 
jectory (see Fig. 1) issues from the point A (0, vo), following the circle 


P+ (xt 


through the point B(r — 1, 0) on the x-axis to D(O, —v). At this 
point the trajectory switches over to the circle 


Z? + (x — 1)? = 7’, 


which it follows back to the point A (0, v9), where it switches again to 
the original circular arc, and so on. The motion represented by this 
trajectory is clearly periodic, with quarter-period equal to the angle 
subtended by the arc AB at its center of curvature C’(—1,0). Thus, 
for the period to have the value 27/w, we need merely require that angle 
AC’B equal r/2w, or that 


= tan (12) 


To find the form of the function x,(t), let P(x, z) be the point on the 
upper arc representing the state of the system at time t, 0 < t < r/w. 
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Since angle AC’P = t, angle PC’B = = — t; therefore, since |C’P| = r, 
w 


ros ( :) and 
2w 


x =x,(t) = (7 -1) (13a) 
2w w 


where 


(14) 


Now let Q(x, z) be a point on the lower arc, corresponding to a 
value of in the interval r/w <t < 2r/w. Since the point D(0, 
corresponds to t = and angle DCC’ = angle DCQ = t — r/w 


and 
angle OCC’ = (: —- ) 
2w w 


Therefore the x-coordinate of Q is given by 


3 
= = | reos (3% - 1) (: <t< *). (130) 
2w w w 


Equations 13, which represent x,(4) for 0 <i < 2r/w, must be 
slightly modified to yield analytic expressions for x,(¢) for all ¢ >0. 
Because of the periodicity, the point A (0, v9) represents the state of the 
system at times ¢ = 0, 27/w, 4r/w, ---, and B(O, —vo) at t = 
32x/w, 5x/w, ---;similarly, the upper arc in Fig. 1 is the trajectory during 
intervals 2kr/w <t < (2k + 1)r/w, k = 0, 1, 2, 3, --- and the lower 
arc during intervals (2k + 1)r/w <t < (2k +2)x/w. Thus, we need 
merely increase by 2kx/w the arguments of the cosine in Eqs. 13, to 
obtain the following representations of x,(t) valid for all t 2 0: 
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2. Restrictions on Parameters 
The preceding construction of x,(¢) requires w > 1, for the angle r/w 


subtended by the arc ABD at C’(-—1,0) is always less than x. If 
w < 1, solutions x,(t) may be found by assuming x, and x, 180° out of 
phase, but we restrict our attention here to harmonic solutions of 
Eq. 4 with w > 1, principally because the mode of construction and the 
form of the results may be conveniently modified for treatment of the 
subharmonic case. 

A further restriction is needed on the parameters yo and w. It has 
been assumed above that x,(t) = x;/(t) + x(t) > 0 for 0 <t < r/o, 
that is, 


2 
x;(t) = sin wt > — x,(t) (0 =). (16) 


Since x, is positive in this interval, the simplest sufficient condition for 
Eq. 16 is that yo be negative, in which case, since w > 1, x,;(¢) will be 
positive. Then x,(t) and x(t) are in phase, but the forcing function 
yow* sin wt is 180° out of phase with them. It is a curious but well- 
known phenomenon in linear systems that the forced oscillation is 
respectively in phase, or 180° out of phase, with the applied force 
according as the natural frequency is greater or less than the frequency 
of the applied force. 

In case Yo is positive, x, will be negative during the first half-period, 
and then appropriate conditions must be imposed on the parameters to 
insure x, = x, +x, > 0 during this time interval. In Appendix B, 
for example, it is shown that, for 1 < w < 2, x,(t) > — x,(t) when 
0 <t < m/w, provided yo is less than some positive quantity depending 
on w. Possible improvements of the result are also mentioned in this 
Appendix. 

The principal results of the preceding analysis are contained in 


Theorem 1. Suppose w > 1. Then if yo < 0, the function x,(t) 
= x,(t) + x,(t) is a harmonic solution of Eq. 4 (that is, one with 
period 2/w), where x,(t) and x,(t) are given by Eqs. 8 and 13 (or 
Eq. 15), respectively. If yo > 0, x(t) = x,/(t) + x,(t) is a harmonic 
solution provided the values of w(> 1) and yo are such that x,(t) 
+x,(t) > Ofor0 <t < w/w. 


The theorem refers to a periodic solution x, which is positive during 
the first half-period. This assumption about the sign of x, was made 
only for convenience ; it is clear that if x, and x, are multiplied by —1, 
xX, = x, + x, will again be a harmonic solution, which is, however, 
negative during the first half-period. 
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Il. SUBHARMONIC CASE 
1. Formal Subharmonic Solution 


We are now interested in periodic solutions of Eq. 4, of least period 
m = 2,3, ---. As in the harmonic case, we prove that such 
solutions exist by constructing them. 

Let m be any positive integer greater than 1. Again we seek a 
solution of the form x,,(t) + x,;(t), where x;(t) is given by Eq. 8 and 
x»(t) is a solution of Eq. 9 with period 2m7/w (obtained in exactly the 
same way as was x,(t) in the harmonic case). Since x,;(t) also has the 
period 2mz/w (it has all periods 2n7/w, n = 1, 2, 3, ---), the function 
Xem = Xm + x, will have the required period. (The subscript sm indi- 
cates ‘subharmonic of order m.’’) 

Since for w > 1, Eqs. 13 furnish a solution of # + x + sgnx = 0 
of period 22/w, a solution of period 2mz/w is obtained merely by re- 
placing w by w/m in Eq. 13. Thus 


1 + sec — cos (= :) (0 <t< 
2w 2w 
Xn(t) = (17) 


mr 3mr mr 
1 — sec cos ( _ :) ( <t<— ) 
2w 2w w w 


when w/m > 1, or w > m. The specification of the function for 
all ¢ 2 0 would be given by equations analogous to Eq. 15. Clearly 
X» Vanishes at t = kmr/w, k = 0, 1, 2, ---, and it is positive for 
0 <t< mr/w, negative for <t < 2mm/w, and so on. It is 
worth noting that the initial conditions giving rise to solution (17) are 
Xm(O) = 0, = vo, where 


(18) 


2. Restrictions on Parameters 


We now know that x,, = x, +x, has the required periodicity 
property, also that both x,, and x, vanish at t = 0, mr/w, 2mm/w, ---, 
but we must make sure, as in the harmonic case, that x,, = X» + X, 
actually satisfies the differential equation (4). In other words, the 
parameters must be subjected to conditions which will insure that x,,, 
is positive during the first half-period, 0 < t < mm/w, negative during 
the second, etc. The question is more complicated here than in the har- 
monic case because x;(t) has (m — 1) zeroes in the interval 0<t<m7/a, 
and not only the values of x; and x,, must be compared, but also the 
values of z,; and 

It will be seen that the conditions to be imposed depend on the sign 
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of yp and the parity of m. The analysis may be conveniently divided 
into two cases: (1) yo < 0, m odd, which requires the weakest assump- 
tions; (2) all the other possibilities. Suppose first that yp < 0 and m is 
odd. In this case both x, and x, start off positive from ¢ = 0 and 
always agree in sign for values of ¢ in the neighborhood of t = kmm/w, 
k =0,1,2,-+--+. As indicated in Fig. 2, x, + x, is certainly positive 


| 


| 


Fic. '2. Comparison of xm(t) and x;(t) for m odd. Solid and dashed curves in upper figure 
represent x/(t) in the cases yo < 0 and yo > 0, respectively. 


for 0<t<2/wand (m—1)2r/w<t<mr/w; furthermore, if x, +x, <0 any- 
where in the interval r/w < t < (m — 1)x/w, it must be in the sub- 
interval < t < 37/2w, where x, takes on its smallest values. 
Therefore x.» = Xm + x; > 0 throughout 0 < t < m/w provided that 
X»(9/w) 2 max |x,|, or, referring to Eqs. 8 and 17, that 


(19) 


Condition 19, then (together with w > m), is sufficient to insure the 
validity of the constructed subharmonic solution x,,, = x» + x, in this 
case (vo < 0, m odd). 

We now turn to the other possibilities (yo < 0, m even; yo > 0, 
m odd or even). For yo > 0 and m odd, x,;(t) is represented by the 
dashed curve in Fig. 2; we see in this case that x, and x,, disagree in 
sign both for ¢ close tot = 0 and for t close tot = mx/w. In Fig. 3 the 
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graphs of x,(t) and x,,(¢) are compared for the case yo > 0, m even; 
the functions agree in sign for values of ¢ close to t = mx/w, but dis- 
agree for t close tot = 0. It is not hard to see that for yo < 0 and m 
even the situation is reversed. 

These three cases may be considered together because essentially 
they all present the same difficulty, which does not occur in the case 
yo < 0, m odd; namely, in the neighborhood of either ¢ = 0 ort = mr/w 
or both, x; and x,, are of opposite sign. To be specific, we shall examine 
| the case yo > 0, m even; it will be evident that the analysis and resulting 
sufficient condition apply as well to all the other cases. Examination 
of Fig. 3 reveals that the critical subinterval is 0 <¢ < r/2w. If we 


Fic. 3. Comparison of x,(t) and x;(t) in the case yo > 0, m even. 


can insure x; +x, > 0, or x, > — x,, on this subinterval, the in- 
equality will hold throughout 0 < t < m/w. 

From the nature of the functions x;(t) and x,(¢), we know that 
—z, and z,, are both monotonic decreasing for 0 < t < r/2w. There- 
fore, it is sufficient to specify that z,,(#/2w), the smallest value of z,, for 
0 <¢t < w/2w, be greater than or equal to —z,(0), the largest value of 
—z, in this interval; then z,,(t) > — #,(t) for 0 < t < w/2w, and since 
Xn(0) = x,(0) = 0, integration from 0 to ¢ yields x,,.(t) > — x,(t) for 
Q <t < w/2w, the desired result. Using Eqs. 8 and 17, we may write 
the sufficient condition in this case in the form 


(m — 1)r 
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Since w > m > 1 and here yo > 0, this is a relation between positive 
quantities. In each case we must compare x,, and |x,| ; it follows then 
that Eq. 20 will apply to all three cases if we write it in the form 
(m — 1) 

i Yo 

(21) 
co 
2w 
It is not hard, indeed, to see that this sufficient condition applies to all 
four cases: Yo positive or negative, m odd or even. In the most favor- 
able case (yo < 0, m odd) we may also use the alternative condition (19). 
We have proved 


Theorem 2. If the positive integer m(>1) and the numbers w and 
yo are such that w > m and inequality (21) holds, then the function 
Xem(l) = x,(t) + x,(t) is a subharmonic solution of the differential 
equation (4) with period 2mm/w, where x,(t) and x, (t) are given by 
Eqs. 8 and 17, respectively. In case yo is negative and m is odd, an 
alternative to sufficient condition (21) is (19). 


The graph of the subharmonic solution x,,,(¢) for the case yp = — 1, 
m = 3,w = 4 is given in Fig. 4. 
x 


2.0 


Fic. 4. Graphs of x;, x.. and subharmonic solution x, = x; + x» as functions of 
twhen yo = 1,0 = 4,m = 3 
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Ill. CONCLUDING REMARKS 
1. Damping 
We have dealt here with a simplified mathematical model of a relay 
servomechanism ; aside from the fact that the on-off element has been 
assumed to be ideal, damping has been neglected. It is intended, in a 
future investigation, to take damping into account. 


2. Stability 


For a given input function y(t) = yo sin wt, Theorems 1 and 2 state 
that if w lies between the integers M and M + 1 there exists a multi- 
plicity of periodic solutions of Eq. 8, with the respective periods 2mz/w, 
m = 1,2, ---, M, provided that in each case yo satisfies an appropriate 
condition. Which of these steady-state solutions arises in a specific 
case (assuming that our undamped servomechanism can really exist) 
depends on the initial conditions x(0), #(0)—more precisely, on the 
value of 


tan 
an =—. 
2w 


These values of v» (for m = 1,2,-+--, M) comprise, however, a 
finite set of numbers, and in practice v) would almost certainly have none 
of these values, thus giving rise to a nonperiodic solution. On the 
other hand, it is not unreasonable to suppose that a nonperiodic solution 
will approach a “‘neighboring”’ periodic solution—a stable one, of course. 
Which steady-state solutions are stable, and whether indeed nonperiodic 
solutions approach them in time, are questions which remain to be 
investigated. In analog computations reported by Ogata (3), however, 
harmonic and odd-order subharmonic oscillations have been observed 
in the case of nonlinearities like the one studied here. 


3. Harmonic Response Relation and Resonance 


In the harmonic case the maximum value of x, = x, + x; occurs at 
t= 2/2w. Thus the amplitude of oscillation, A, = max |x,|, is 


given by 
us 
( 2w 


Suppose, to simplify matters, that yo <0. Clearly both terms on the 
right get unbounded as w—1. We have, then, the possibility of 
resonance, as in the case of a forced undamped linear oscillator. 

When w» gets large, however, A, — — yo, that is, the response is 
purely a “‘forced’’ one, whereas the response of an undamped linear 
system always contains forced and free components. 
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4. Subharmonic Response to High Frequencies 


Suppose w is much larger than 1, and | yo! is small. Then max|x;| 
is small, and, as we have just seen, the harmonic response amplitude 
A, isalso small. It is possible, however, in this situation for the system 
to exhibit a subharmonic response of order m (m = 1, 2, ---, M < a). 
If mis rather large, max |x,,| is not negligible ; then the amplitude, A,n, 
of this subharmonic response ¥,,, = X, + X; is approximately given by 


Aim = max |x,,| = sec — 1. (w >m> 1) 


As m/w — 1, A,» gets unbounded. This possibility of exaggerated, and 
perhaps destructive, response points up the practical importance of 
analyzing the various types of oscillations of which a nonlinear system 
is capable. 


REFERENCES 


(1) R. L. CosGrirr, “Nonlinear Control Systems,’ New York, McGraw-Hill Book Co., Ine., 
1958. 

(2) J. C. West ann J. L. Douce, “The Mechanism of Sub-Harmonic Generation in a Feed- 
back System," Proc. 1.E.E., Vol. 102, Pt. B, p. 569 (1955). 

(3) K. OGata, “Subharmonic Oscillations of Nonlinear Feedback Control Systems," Paper 

No. 58-IRD-7, Am. Soc. Mech. Engineers. 


APPENDIX A 


Consider a simple servomechanism with proportional plus discontinuous control, whose 


governing differential equation is 


- Tsgne 


definition 


By 


(A-2) 


sgne =signume = 


(e < 0) 


Let us interpret Eq. A-1 in terms of a specific physical example: a mechanical system in which 
the input and output are rotations of shafts. Then the symbols in Eq. A-1, and their units of 


measurement, may be defined as follows: 


r = time, in seconds; 


r = output (angular displacement of a shaft), in radians; 


¢ = input (angular displacement of a shaft), in radians; 


¢ =r —c¢ = error, or deviation of output from input, in radians; 


I = constant (moment of inertia) of output element, in slug-ft?; 


k? = constant (torque per unit angular displacement) of proportional control element, in 


lb-ft per radian ; 


constant (torque) of discontinuous control element in lb-ft. 


3 
Bere 
d* 
dr 
+1 (e > 0) 
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The constants k*, J and T are, of course, positive. 
If we set r = ¢ + cin Eq. A-1, we obtain 
de d* 


of. 2p 4+. 7 ao = 
+ ke + Tsgne (A-3) 


By introducing the frequency parameter 


and the dimensionless quantities 


i = vr, 


we may write Eq. A-3 in the simpler form 


F+xu+sgnx = — j, (A-6) 


where - = d/dt and -- = d*/df@. Sgn x = sgn e, that is, x and e are of the same sign, because 
k?/T is positive. The variables /, x, y and s are referred to in the paper as the time, error, 


input and output, respectively 


APPENDIX B 


We wish to find a condition on » and ¢ which insures x,(1) > —x,(1) for 0<t<r/w(w>1), 


where 


= 1 + cost + tan sint, — x/(t) = c sin wt 


For convenience we have set 


1 w 


Note that because w > 1, yp and ¢ have the same sign. 
First of all, since the graphs of x;(¢) and x,(t) are symmetric about the line ¢ = x/2w, we 
need only deal with the interval 0 <¢ < #/2w. From examination of the graphs of sin @ and 


cos @ it is readily seen that 


2 — 2 
sing <¢@ and cos > - 1 


for 0 <@ < 80 that for0 <t< r/lw 


2 x\2 2 2 © 
xr, > -1l--t4+14 tan (5 ) —-t+-—ftan ; 
2w 


and 


It will be true, then, that x, > — x; for 0 <t < r/lwif 


2 
cut < ( -1+ tan 

2w 

2 
cs (- 1 + tan 

yar) 


For w < 2 this upper bound on ¢ (which corresponds to an upper bound on yo) is positive. 
For w > 2 it is negative, but we know that we can always do better than this, for any negative 


or 


2 


c (or yo) will do (see Part Il, Sec. 3) 


k 
(A-4) 
v/ 
7 
cz 
— xy < Cul, 
= 
| 
: 
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It is possible to improve this result by using sharper estimates of sin @ and cos ¢, ob- 
For example, if we integrate 


tained by integrating the inequalities given above. 
2 

@e<sing<@ (0 << 
< 


from 0 to @, we get 


<cos@ <1 (0 <#<3) 
2 


¢ , respectively, and it is possible to show, 


The amplitudes of x; and x; are r — 1 and 
indeed, that there is a value of w, wo, such that for all # > wo > 1, the condition ¢ <r — 1 


implies x,(t) > — for 0 <t < w/w. 
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CONTRIBUTION TO THE MODIFIED Z-TRANSFORM THEORY * 


BY 


E. I. JURY! 


SUMMARY 


In this discussion, a useful theorem applied to the modified z-transform (1, 2)? 
is introduced and verified. This theorem is fairly general and covers a few special 
cases that found many applications in the study of pulsed-systems. 

These mentioned cases are discussed extensively and documented with examples 
and a short table in the Appendix. Various other cases of the theorem and its ex- 


tension are pointed out in this paper. 
CONVOLUTION THEOREMS 


Theorem 1. Complex Convolution Theorem for the Modified Z-Transform'® 


S(t) X h(t) = p-'F*(p, m)H* m) dp, 


O<m <1, 


I 
m) F* m) dp 


f(t) = h(t) =0 for t <0 


where IT is a closed contour in a positive sense in the p-plane which 
encloses all the poles of p-'F*(p,m), and F*(p,m) = jn{ f(t) J--», 
H*(z,m) = h(t) 

To prove the above theorem, we apply the complex convolution 
theorem to obtain the Laplace transform of multiplication of two time 
functions as follows (1-4) : 


T 
LU f*(t, m) X h*(t, = f F*(q, m)II*(s — q,m)dq (2) 
7) 


where the line integral ¢ — j r° toc +j r° lies on an analytic strip in 

* This research was supported by the United States Air Force through the Air Force Office 
of Scientific Research of the Air Research and Development Command, under Contract No. 
AF 18(600)-1521. Reproduction in whole or in part is permitted for any purpose of the United 
States Government. 

' Department of Electrical Engineering, University of California, Berkeley, Calif. 

* The boldface numbers in parentheses refer to the references appended to this paper. 

* The complex convolution theorem for the z-transform case has been developed by 


Cypkin (2). 
114 
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the complex g-plane that does not pass through or enclose the poles of 
F*(q, m) or of IT*(s — m) and 


F*(q,m) = f*(t, (3) 
H*(s,m) = £[h*(t, m) }. 
The typical function f*(t, m) is defined as 


f*(t,m) = — (1 — m)T] X Sr(t) 
= f(nT — T + mT)é(t — nT) O<m 


n=O 


The Laplace transform of f*(t, m) and h*(t, m) is our definition of 
the modified z-transform of f(t) and h(t), respectively, thus 


F*(s,m) = £[f*(t, = S(O], 
LLA*(t, m)] = jmLh(t)] = H*(s,m). (6) 


Because of the periodic nature of 67(t), the Laplace transform 
F*(s, m) is a rational function of e™; thus, writing F*(s, m) in terms 
of replacing e7* by z we denote this function as F*(z, m) | 

Changing the path of integration in Eq. 2 from the g-plane into the 


de’ 


ere 
p-plane (1), (where p = e7¢), and writing for dg = Tet’ we obtain: 


LL f*(t, m)h*(t, = p-'F*(p, m)HT* m) dp. (7) 


It should be noted that the line-integral in Eq. 2 is changed into a 
closed contour T in the p-plane as required by the change of the variable 
of integration. From the definition of the modified z-transform, the 
left side of Eq. 7 is equal to the 3n[_f()A(t) )]. Therefore : 


f(t) X = sade? 'F*(p, m)H* m) dp 


which proves the theorem. 


For the case h(t) = f(t), 


= p-'F*(p, m) F* m) dp. 


‘When A(t) = y(t), the left side of this equation reduces to F*(2, m) 


| 
| 
| 
=, 
| 
| 
| 
cr: 
| 
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Evaluation of the Complex Convolution Integral (1, 2) 

To evaluate Eq. 8, we can make use of Cauchy's Integral Formula 
and evaluate the residues inside the circle whieh encloses only the poles 
of p-'F*(p,m). An alternate procedure is to obtain the residues of 


H* (2 m) outside this circle, with path of integration taken in the 


negative sense, provided the integral isconvergent. The choice between 
the two methods depends on the number of poles of each of the functions 
F* and 

Making use of the first procedure and assuming /*(p, m) has simple 
poles, then 


A(p,, m) ( 
mL S(t h(t)} = H* 10)* 


where 
A(p, m) 
B(p) 
B(p) Pi, Pe, 
dB 
dP | pay, 


In general for multiple poles cases, Eq. 8 becomes : 


X J] = Residue of p-'F*(p, m)H* m) 


at poles of 


p 'F*(p, m) 


(11) 


F* (p,m) = 


(12) 
B’(p,) 


ILLUSTRATIVE EXAMPLES 
Example 1 


Assume 
f(t) = e-*'| 
h(t) = e-#*| (14) 


The modified z-transform of each of the above functions can be ob- 
tained by using the following relation or through the use of tables (1, 2). 


F*(z,m) = ef f[t — (1—m)T] xX 
= > f(n—1+m)Ts-*, O<m<i1 (15) 


Thus from tables (1) or from Eq. 15, 
e amT 
F*(z,m) = (16) 
5 For these cases it is assumed that the residue of the pole at the origin is zero. If such a 
residue exists, then it should be modified to include its value. 


4 
: 
4 
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H*(z,m) = = 
From Eq. 8 


= 


xj rp 


Since F*(p, m) has a simple pole at p = e-*7, Eq. 10 can be applied 
which yields: 
e~ (a+B)mT 


e-aT — -8T 


Example 2 
Find the modified z-transform of the following function (2, 5) 
(t) 


0 
1)? (2 ) 


f(t) ] = F*(z, m). (21) 


In this case (1): 


From Eq. 8 


Tp — m1 
f(t)] = “(2 m) dp. (22) 


The function of the integrand in the region of integration has a 
multiple* at p = 1, thus the sum of the residues equals (1), 


| 


4. 


(23) 


* Since in most cases, (5, ml = (), it is assumed that the residue of the pole at the 


origin is zero. 
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of 
It can be easily verified that for all forms of F* (Z m)| , the 


following is true 


(24) 


Therefore, 


0 
ef (t) = r| (m — 1)F*(z,m) F*(z, m)| = F,*(z,m). (25) 


Similarly for 
a 
f(t) = 7 | — 1)F\*(z,m) — F,*(z, m)| (26) 


By iterative extension of the above one can obtain, 4,{t* f(t) J, k = posi- 
tive integer. 
APPLICATIONS OF THE THEOREM 


Several special cases of wide interest can be obtained from the men- 
tioned convolution theorem : 


(a). Application to obtain 3f f(t)h(t) ] 


The z-transform can be obtained from the modified z-transform by 
multiplying by z and letting m = 0 in F,*(z, m). 


F,*(z) = sF,*(z, m)| moo = }| mao. (27) 


Alternately, in integral form /,*(z) can be represented (2): 
) 


* = = * 
F,*(s) = 3Lf(th(t)) > x H > dp 


F*(p) = 3L f(t) 
H*(z) = 3[h(t)]. 


where 


2 


(b). Application to obtain [f(nT)}? 


n =() 


Letting h(t) = f(t) in Eq. 28, we obtain 


n=O) 


Op p p=! 02 
; 
‘ 
(28) 
(29) 
(30) 
(3 | ) 
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Letting z = 1 in the above equation to obtain 


Lf(nT)} = F*(p) F*(p-')p-'dp, (32) 
2rjJr 

where the contour-I° is along the unit circle in the p-plane, because all 
poles of p~'!F*(p) should be inside the unit circle in order for the summa- 
tion to exist. 

The above equation is used extensively in evaluating the sum of the 
square error at the sampling instants in sampled-data control systems 
(2) and also in obtaining the noise power gain (6). The evaluation of 
the above integral (without obtaining the poles of F*(p)) can be ob- 
tained as ratio of two determinants which has been tabulated in some 
publications (2, 5). 


(c). Application to obtain (5). 


f(HA(t)dt = F*(p, m)IT*(p-', m)p~'dpdm (33) 
£7) 0 


f [ f(t) = J, F*(p, m) F*(p-', (34) 


where the closed contour I is the unit circle in the p-plane. 

The above theorem is analogous to Parseval’s theorem for Fourier 
transforms. It has been applied for the analysis of integral-square 
error in sampled-data control systems (5). 

To show the above application from the general theorem we first 
obtain the following equivalent form of the integral in the left side 
of Eq. 33. 


n=9 79 


= (n—1+m)T, n = integer O<m<1 


On the other hand, the left side of Eq. 8 can be written as 


f(t)h(t) = f(nT — T + mT)h(nT — T + mT)z-". 


nomi) 
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If we let z = 1, multiply by 7 and integrate w.r. to m we get: 


rf 
= f(nT — T + mT)h(nT — T + mT)dm. (37) 


Interchanging the order of summation and integration and using Eq. 8, 
we obtain 


m)H*(p-', m)p—'dpdm. (38) 
2xj 


Ty 


JO 


From Eq. 35, the above equals 


0 


i 
f(t)h(t)dt = fro. m)H*(p—, m)p-'dpdm. (39)? 
r 


When h(t) = f(t), the above yields Eq. 34. The evaluation of this 
integral is simplified if the order of integration in z and m is inter- 
changed.* This is possible since the integrand is continuous along 
and in the interval 0 < _m <1. The evaluation of the first integral (5) 
is carried similar to Eq. 32. 

Noting Example 2, for ;,.[¢/(t) ] and letting h(t) in Eq. 39 be [tf(t)], 
we get 


1 
f f(t) Pdt = — 1)F*(p, m) F*(p, m)| 


x F*(p-', m)p-'dpdm. (40) 


Similar relations can be obtained for higher time moments. 


(d). Application to obtain 


T p-'F*(p, m)H* (: m) dpdm 
— 1) Jo Jr p’ 
F*(p, m)H m) dpdm. (41) 


7 Asa "spec ial case when h(t) = p(t) = unit step, we obtain 


f = rf F*( p, m) dm, where F*(p,m) = f(t) 
0 p=! 


* See Theorem 2 for evaluation of the integral w.r. to m. 
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The above relationship can be obtained from writing the integral in a 
summation form as follows: 


= r|  f(n — 1+m)Th(n — 1+ m)Tdm 


J 


0 


+ — 1+ m)Th(n — 1+ m)Tdm|. (42) 


The modified z-transform of the above yields® 


f(n — 14+ m)Th(n — 14+ m)Ts-"dm 


+T fin— 1+ (43) 


n=) 


but the infinite summations are the modified z-transform of f(t) A(t), 
thus: 


= - F,*(z, m)dm + rf" F\* (z,m)dm (44) 


where 


F,*(z,m) = f(HA(t)) = f F*(p, m)p-'H* m) dp. (45) 


For the special case when h(t) = w(t) = unit step, the above yields 


t 
= F*(c,m)dm + rf F*(z,m)dm_ (46) 


where 


F*(z,m) = (47) 


Equation 46 also yields the modified z-transform (2) of ; 
§ 


Exten- 


sion to obtain 3» _ |, by iterative process is straightforward. 
§ 


(e) Application to obtain 


F(s) = T m)dm, 


Noting that (1), F%(2) = 


| 
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where 
F*(z,m) = jn{f(t),t F(s) = £[f(t)). (48) 


To show the above, we multiply the left side of Eq. 33 by e-*! and inte- 
grate to obtain 


if 


— T + mT)h(nT — T + 


(n—1+m)T, O<m<1 (49) 


or 


f(t)h(the-*'dt 


= met f(nT — + mT )h(nT — T+ mT)2-"dm (50) 


n=( 


but 


x f(nT — T + mT)h(nT — T + mT )2-" = jal f 


therefore from Eq. 8 we get 


F,(s) = f(h(tHe-*'dt 


0 I 


1 
Fi(s) = rf m)dm, 


where 


F\*(2,m) = jnLf (h(t) and Fy(s) = (53) 


The above relation shows the equivalence between the Laplace 
transform and the modified z-transforms. ‘The above has been shown 
derived differently by Cypkin (2).. It is obtained here in a simple 
manner as a special case of the convolution theorem. 


t F*(s, m) can be directly obtained from F(s) by using the following convolution in- 


tegral (1): 


| 
F*(s,m) = F(p) — dp 


] 
|| 
= 
4 
n=O 
(52) 
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Example 


amT 
Knowing F*(z,m) = , obtain F(s). From Eq. 52, Fis 


equals 


F(s) 


Ts f oe rat 
= mr = 5 
Z—e-*T Jo s+a (54) 


~ 


(f) Application to obtain 


e~j(r/T) 


s) == ‘here = 
F(s) F (q) where F*(q) = F*(s)|,24. (55) 


Knowing F*(s) = £[f*(t)] from the above equation, to obtain a 
continuous function f(t), whose values coincide with f*(t) at the sam- 
pling instants only, the following condition should be satisfied : 


= u(t) = unit step (56) 
or in Laplace transforms, 


F(s) = £[f(t)u(t)) = 


Using Eq. 2 to evaluate the right side of Eq. 57, we obtain (2) 


e+ j(r/T) 
F(s) = (58) 
s—q 


e~jlr/T) 


In evaluating the above equation only the poles in the linearized strip 


between ¢ — 7 


.toc +] 7 are considered. 


7 


Example 


f*(t) = e-*"T5r(t), F*(q) = F*(s) = (59) 


— 


Evaluating Eq. 58, for the poles of F*(g) we obtain 


F(s) = T— — = 7 (60) 
s+a 

eT] 

dq 


: 
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Alternately evaluating for the pole g = s, we obtain 


1 1 


(g) Application to obtain 


1 
F*(z, m) = 2 F*(2,e27 m), Zn = e(T (62) 


n—l 


F*(z) = > (63) 


The above theorem signifies that, knowing the modified z-transform 
with respect to the period (7/n) (when n integer is larger than unity), 
one can obtain from it the modified z-transform of F*(z, m) with respect 
to the period 7. This derivation is of importance in the analysis of 
multi-rate sampled-data feedback systems (1, 2). 

The inverse of Eq. 62 indicates a time function which is formed by 
the product of the inverse of F*(z,, m) and a sampled step-function with 
period 7, that is, f(k, m)T = f(k,m)T/n X w(kRT). Thus using the 
convolution integral we obtain: 


1 
on 
F (z, m) (Pu, m)p, 1 (64) 


where 
F*(p,, m) = F*(z,, m)|, <5, = £Lf*(t, = S(O) (65) 


and 


= u(t) 


Equation 64 can be written as: 


1 
F*(z, m) F*(p,, m)p, ip. (67) 


)” 


Integrating in a negative sense w.r. to the poles of TMS. jp Which are 


= k = 0,1,2, ---,2 — 1, 


(J. F. 1. 
: 
5 
F(s) = = (61) 
: 
or 
: 
: 
1 1 
= . (66) 
1 — | 1 — we Tine =3." 
: 
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we obtain the sum of the residues 


=— F*(p,, m)p,~' d . (68) 
dp (p./8.)* — 1 "se 
Thus, 
F*(z, m) = m), (69) 


k= 


For the case of z-transform, letting m = 1 in Eq. 69, we obtain 


F*(z) = (70) 
k =f) 


Further application of the theorem in the analysis of modulated 
discrete signals is of further use in the study of modulated sampled-data 
systems. 


ADDITIONAL THEOREMS 


Theorem Ja 


lim f(t) = 0, t=(n—-1+m)T. (71) 


This theorem can also be derived from Eq. 8 as a special case when 


h(t) = r However in using the convolution integral it is necessary to 


know jm [? | To avoid this, the above is derived differently as follows : 


F*(s,m) = f(nT — T + = f(nT — T + mT)s-*, (72) 


n=) nel 


since it is assumed that 


S(t) | = 9. 


Multiply the above equation by z~” and integrate to obtain 


m)dz = f(nT — T + 


—-n—-m+1 


O<m<1 (73) 


: 
5 
| 
i 
. 
d 
sie 
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(nT —T+mT) 
1 + m)T 


If f - 2” is added to both sides of Eq. 74, then this term 
(n —-1+m)7 
has a contribution only at m 


f(t) 
= 1, which is the lim ——, thus the right 
(t 
side of the above equation is 3, f | 


.| = 72 F*(z, 


f(t 
) <m < i. 


“0 
For the case of z-transform we obtain, by substituting m 


= 1 in Eq. 
| F*(z)dz + 


n 


where 


F*(z) F*(z, m) | 


Extension of Eq. 75 to obtain .| 


| is evident (2). 
Special cases which stem from this theorem are indicated in the 
following. 


(a) f(t) | 
2 == 
"Li+aT 


where / is a positive constant larger than zero. 


F*(s)dz 


The above is obtained 
from Eq. 75 by substituting for m — 1 = h, and using the z-transform 
case. 


(b) When z 


= 1, relation (a) yields: 


x 


f(nT) 1 
(+h) F*(2)dz, h> 0. 
(2 + h) T 
For the case h = 0, the above is modified to yield 


x 


f(nT) 1 


z-'F*(z)dz + lim f(nT) 


(nT) 
When m = 0, we obtain ) = 


or, 
: 
: 
(75) 
(76) 
> 
~ ~ ( ) 
2 
; 
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Example 
Obtain, using Eq. 75 


F°(s,m) = — 


~ 


me—amT 
| as 


for0 <m <1 


z-"F*(z, m)dz 


Finally, 


1 — eer 2 


bd aie -2aTp-amT ( < < &3 
+ in »+ +, (83) 


For the z-transform let m = 1 in Eq. 83, to obtain after closing the series 
(84) 


A short table of 2- transforms i is give en in the Appendix. 


"It is of interest to note that we 


For the case m = 1, we obtain 
= &{f) + ain = lim | ]. 
)+ where (t im [ ; | 
In some publications (2) the z-transform of [7] is defined only for ¢ > 0, thus the infinite 


value of 6,(t) is ignored. 


5 
127 
: 
t s-¢ 
2-"F*(z, m)dz = 
(81 
1 2 
m m + m + 
m mati m+2 
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Theorem 2 
rf F*(z, m)H*(z-', m)dm = 3(F(s)H(— s)], O<m<1. (85)" 


The above theorem is an extension of the following theorem which is 
discussed in the literature (1, 5, 7, 8). 


F*(z, m)F*(z—', m)dm = 3[ F(s) X F(— s)]. (86) 


Following similar steps for the derivation (5) of Eq. 86, one can 
easily obtain the theorem in Eq. 85. 

The use of the above theorem is of advantage in the evaluation of the 
integral w.r. to min Eq. 39. Thus from Eqs. 39 and 85 we obtain 


f = ff F*(p, m)H*(p—, m)dmp-'dp 
| 


CONCLUSION 


The complex convolution integral for both the z-transform and 
modified z-transform has been discussed in detail. The applications of 
the theorems for various fields such as sampled data systems, digital 
control systems, circuit theory and to the operational solution of differ- 
ence equations both linear and non-linear, are possible using the tech- 
niques of this paper. It is of particular importance that the extension 
of the convolution integral for solving certain types of non-linear 
difference equations is warranted. 
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is [F(s)H(— s)]. 
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APPENDIX: A SHORT TABLE OF Z-TRANSFORMS OF CERTAIN FUNCTIONS 


F*(s, m) F*(z) 
aT 
at U! amt") a>0 


+ 


Zz 


DT 


” 


sin at 1 ifr z—cosaT 


2x 
w,=—, k=integer 


T? f 2 F°(s, m)ds+lim 1 \ds+lim ; 
O<m< 

F*(z) 


T 1 


k>0, and F,*(s, m) = F,*(z) = 3[t*- (1) 
integer 


: 
fit) 
129 
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THE FIRST ELECTRIC COMPUTER, A MAGNETOLOGICAL 
ANALYSIS * 


BY 


GEORGE W. PATTERSON ! 


SUMMARY 


Mechanisms designed to carry out the ordinary operations of decimal arithmetic 
were invented several centuries ago, and indeed Charles Babbage, in England, con- 


ceived automatic sequence control over one hundred years ago. The modern elec- 
tronic high-speed computer based on binary Boolean logic has other ancestors besides 
decimal mechanical calculators in its pedigree, however. The mechanical realization 
of binary logic is generally credited to William Stanley Jevons, British economist and 
logician, in 1869.2. This mechanism is an ancestor of the modern computer on the 
binary side of the family. It has recently come to light that an American university 
professor designed what is probably the first electric binary logic machine, and this 
is very likely the first design of an electric data processor of any kind 


Professor Alonzo Church of Princeton University has discovered a 
cardboard schematic diagram among the papers and manuscripts of 
Allan Marquand (1853-1924), former Professor of Archeology and the 
History of Arts at Princeton University. This diagram, believed to 
have been prepared about 1885, was described and reproduced in 
facsimile in a communication to the magazine Science... Marquand 
had constructed an earlier mechanical logic machine‘ that still is extant. 
It seems quite plausible that the electrical design was based on the 
earlier mechanical one, but it is not clear whether the electrical machine 
was ever built, or, if built, whether it operated successfully. The 
present writer would give the designer the benefit of the doubt, and 
believes that the proper values of the parameters (not given on the 
schematic), if not known initially, would be soon discovered. 

The explanation of the mode of operation of the machine is not 
straightforward since a number of assumptions must be made concern- 
ing the magnetomotive forces operating on the relay armatures. The 
basic difficulties arise from two facts: (a) The truth-functional opera- 
tions are not mechanized in terms of controlled-contact networks, but 


* This investigation was performed under the sponsorship of the Department of the Army, 
U.S. Army Signal Research and Development Laboratories, Contract No. DA36-039 SC-72344. 

' The Moore School of Electrical Engineering, The University of Pennsylvania, Phila- 
delphia, Pa. 

? MARTIN GARDNER, “Logic Machines and Diagrams,’ New York, McGraw-Hill Book Co., 
Inc., 1958, p. 91. 

*W. Mays, “The First Circuit for an Electrical Logic Machine,"’ Science, Vol. 118, pp. 
281-282 (1953). 

‘ALLAN MARQUAND, “A New Logical Machine,"’ Proc. Am. Acad. Arts and Sci., Vol. 21, 
pp. 303-307 (1885-6). See also Gardner, op. cit., pp. 106-113. 
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in terms of threshold elements whose logical properties are described 
in a recent paper*®; (>) the relays are not exclusively truth functional 
in their operation, but also realize storage (memory) devices. These 
two points will be considered in more detail. Concerning (a), the 
description of a threshold element is as follows 


J 
h=t(Xc(fi) + 2 


i=l j=l 


In this formula, 4 names the proposition : 
the output is in its fiducial state, 

while f; and g; correspond, respectively, to propositions of the form: 

the 7** excitatory input is in its fiducial state, 

the j" inhibitory input is in its fiducial state. 


By fiducial state, we mean that a signal is present. The integral con- 
stants, J, J, T correspond, respectively, to the total number of excita- 
tory inputs, the total number of inhibitory inputs, and the threshold. 
The bar or vinculum over the g; serves to negate the proposition, that 
is to say, 9, will be true if and only if g; is false, that is, the 7‘ inhibitory 
input is not in its fiducial state, that is, it is in the antifiducial (no signal) 
state. The expressions c(f,;) and c(g,) denote characteristic functions 
of the propositions. c(f/,) is 1 if f; is true and is Oif f; is false. From 
what has just been said concerning the bar, it follows that c(q;) is 1 if 
g, is false and is 0 if g; is true. 

It turns out then that > c(f,) is simply the total number of excita- 


J 
tory inputs out of J that are in, and ¥ c(g,) is the total number of 


j=l 
inhibitory inputs out of J that are not in, their fiducial states. 

T is the threshold. If the total number of inputs that are either 
excitatory and in, or inhibitory and not in, their fiducial states is greater 
than or equal to 7, then and only then is/A true. For example, suppose 
T is 5, Tis 4 and J is 2. Then if all four of the excitatory inputs are 
energized, and one or both of the inhibitory inputs are not energized 
the output will be energized. 

In the present design, the inputs are physically realized by the five 
windings on the relay core. According to the diagram in Fig. 1, these 
windings all contribute ampere-turns of the same sign. This assumes 
that, and there is no reason to the contrary, the diagram accurately 
reflects coil and battery polarity. In these circumstances it is con- 
venient to assume, for each literal winding, that its fiducial state corre- 
sponds to current flowing in it, and that the fiducial state of the relay 
is the energized state (armature closed against core). Then all the 


+A. W. Burks aNnp I. M. Cort, “The Logical Design of an Idealized General-Purpose 
Computer,” Jour. FRANKLIN INsT., Vol. 261, pp. 299-314 and 421-436 (1956). 
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literal inputs are excitatory, J = 0, and the relay realizes a symmetric 
truth function, one that is true provided 7 or more inputs are energized 
and false otherwise. 

Threshold elements of this simple type can physically realize, 
among other things, generalized n-input ‘“‘or’’ (alternation or disjunc- 
tion, with JT = 1) and generalized n-input “and” (conjunction, with 
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Fic. 1. Diagram of electrical logic machine. 
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T = n), but cannot realize negation. For this reason it is necessary to 
provide the machine with a pair of inputs for each variable. If the 
variable is affirmed, one key is depressed, if denied, the other key is 
depressed. There is no physical bar to depressing both keys simul- 
taneously, leading to a tautologous (or contradictory) term. If neither 
key is depressed the term does not contain the variable. This is an ad- 
vantage since it permits entering a formula in a normal form that is not 
canonical (Boolean, fully developed about all variables, distinguished). 
It is clear that different windings may contribute different amounts 
of magnetomotive force. It seems likely that all windings connected 
to literal keys produce the same number of ampere-turns, say /. We 
use this as a unit of magnetomotive force in the sequel. The fifth 
winding on each core is connected to the operate switch, and is differ- 
ently arranged. The switch produces two distinct non-zero values of 
mmf (contacts 1 and 2 in Fig. 1), as well as zero mmf (contact 3). 
What is the relative effectiveness (O;, O:) of these positions com- 
pared to the other four literal windings? A natural assumption is 
that position 1 is twice as effective as a single literal winding, (O, = 2/) 
since battery B,, which is connected to contact 1, appears to have twice 
as many cells as battery B,, which is connected to the literal keys; for 
simplicity assume that O, = /. These assumptions are believed un- 
workable for reasons yet to be explained. Consider the magnitude of 
the threshold 7 in the light of (b) above: the relays are not exclusively 
truth functional in their operation, but also realize storage (memory) 
devices. The input necessary to cause a relay to close is greater than 
that required to retain it in the closed condition. If T is the operate 
threshold, and the imminent release point is t, we may have 7/t as 
large as 10. Thus if 7 or more is required to cause the relay to pull in, 
we cannot state that less than 7 will cause the relay to drop out. 
Excluded middle is violated and the state of the relay for inputs be- 
tween ¢ and T is indeterminate and depends on the past, and this permits 
data storage or memory. The variable airgap in the magnetic circuit 
produces a double-valued function of flux versus mmf, and more im- 
portant, a double-valued function of armature position versus mmf. 
This is analogous to the hysteresis loop in a rigid magnetic structure and 
the double-valuedness is what makes the memory or storage function 
possible ; the quiescent point is located in the double-valued region. 
The value of the ratio 7/t must be established. According to 
Mays, ‘Position 2 sustains the pointers at the horizontal . . . but does 
not raise the fallen pointers.’’ Also, ‘‘Position 1 raises all the fallen 
pointers to the horizontal.’’ The assumptions are acceptable and yield 
t<¢O,.<T <¢0O,;. If we retain the first assumptions concerning the 
values of 0, and 0,2, we obtain ¢ < 1 < T < 2l, and simple assumptions 


would be to set t = 7 T= 2 to give ample margins of operation. These 
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assumptions yield 7/1 = 3, which is conservative. Unfortunately, no 
method of successful operation has been discovered for this set of as- 
sumptions. Now let us suppose that t < O, = 1 < T < O, as before 
but lett <1 < 51 < 7. This gives 5 < T/t, which is physically still 
feasible. To give wide operating margins we might set ¢ = 0.5, 
T = 5.5 with 7/t = 11, although this is a rather extreme figure. 

The machine is conveniently viewed as a truth-table schema con- 
sisting of sixteen cells, or atoms, upon which we wish to impress an 
arbitrary four-variable truth function expressed in (not necessarily 
distinguished) conjunctive normal form. Each conjunct is the dis- 
junction of eight or fewer literals (the variables, negated and un- 
negated). Usually a variable will not appear both negated and un- 
negated in a disjunction, and as stated before one of the limitations of 
simple threshold elements of the type just described is that they cannot 
realize negation directly. Hence the absence of an A signal does not 
imply the presence of the a (that is, a not-A) signal. This also has 
the advantage that conjuncts may omit variables entirely, and this 
would not be possible if we had four two-position switches for the vari- 
ables, instead of eight one-position switches. 

We successively impress disjunctions of literals on the machine, 
seriatim, yielding a truth table that corresponds to a conjunction of the 
new disjunction with the result obtained from the previous steps. If 
we start with a tautology (corresponding to a conjunction of zero con- 
juncts) then the accumulated result will represent the truth table 
obtained by conjoining all the disjunctions with a tautology. This 
process must be monotone in the sense that each conjunct can change a 
relay indication from energized (truth) to de-energized (falsity) but 
not vice-versa. The explanation is recursive. The starting configura- 
tion is obtained by putting the operate switch on 1. This produces on 
all coils an mmf O, such that 7 < O, and all sixteen relays operate and 
indicate truth. The operate switch is then instantaneously (make- 
before-break contact) placed on 2, yielding an mmf O,. such that t < O2, 
and all relays continue to indicate truth unless O, is removed. This 
manipulation impresses a tautology (corresponding to a conjunction of 
zero conjuncts) on the machine, with all relays indicating truth. 

Now consider an arbitrary truth table on the machine with operate 
switch at position 2 (mmf O,); some of the relays indicate truth and 
some falsity. We depress some (non-void) combination of the literal 
keys. (Note that coils on the A relays are energized by the a key, etc.) 
This applies an additional magnetomotive force of at least / on each 
coil that corresponds to one or more literals. In other words, a force 
of / (or greater) is applied to all atoms that imply the disjunction, and 
to no others. An atom indicating falsity (de-energized relay) will not 
be changed to indicate truth, however, since the mmf is at most 5/ < T 
and the relay remains de-energized. The operate switch is now set to 3. 
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All the atoms that do not imply the disjunction are de-energized and 
indicate falsity. In this way the only atoms remaining energized and 
indicating truth are those corresponding to the conjunction of the 
previous result with the current disjunction. (These disjunctions are 
called premisses by Mays and also by Marquand himself, who views 
them as negations of conjunctions.) It remains only to store the new 
truth table. This is accomplished by returning the operate switch to 
position 2. The magnetomotive force is the same as before but with the 
significant difference that zero or more relays have released. We may 
now release the literal keys, and the storing mmf O, retains the result 
until another disjunction is conjoined. If for any reason we fail to 
retain any literal keys in the depressed state while setting the operate 
switch to zero, all cells go to zero as might be expected since we are con- 
joining a void disjunction, and the latter corresponds to a contradiction. 
A comparison will reveal that this analysis of the behavior of the 
machine is different from that described in Mays’s paper. No engineer- 
ing assumptions consistent with his description and at the same time 
yielding a logically workable machine have been discovered. In order 
to discuss Mays’s description we must translate his vocabulary (con- 
sistent with Marquand’s) intoourown. We use disjunction (a conjunct 
in the conjunctive normal form) for Mays’s premiss. 
Marquand considers the premisses to be denials of conjuncts. He 
expresses each premiss in the form 


o(A, B,C, D) => F, 


where 0 is some truth-functional operator in four or fewer variables 
equivalent to a conjunction of literals. The premiss itself is equivalent 
to ~ o(A, B,C, D). 

In Mays’s analysis of the machine the operate switch is first set to 1, 
and then the individual letters of the first premiss are depressed seriatim. 
These depressions are interlaced with pairs of manipulations of the 
operate switch, first to 3 and then to 2. 

Suppose the first premiss is A = B, which is equivalent to ~ (Ad). 
In our terminology this is the first disjunction, a V B, of a conjunctive 
normal form that may contain other conjuncts. In general, we may 
assume that the machine has an arbitrary truth-functional configura- 
tion already impressed upon it, say 0*(A, B, C, D), and the operate 
switch ison 2. With the operate switch resting on contact 2 those and 
only those relays consistent with o*(A, B, C, D) are operated, although 
all are drawing the same current since no literal keys are depressed. 
According to Mays, the premiss is recorded on the machine in the 
following way. First the A key is depressed. This sends additional 
current thru the a coils (but not enough to pick up any of the relays).® 


* The parenthetical remarks are applicable to the present analysis of the mode of operation 
of the machine. 
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The operate switch is then set to position 3, and only the a coils now 
receive any current. (All non-a relays that had previously been 
operated drop out, leaving only o*(A, B, C, D) & a impressed on the 
machine). The operate switch is then set to position 2. (This restores 
the holding current, storing 0(A, B, C, D) & a, and the A key may now 
be released, since it cannot ever pick up any relays.) Literal key } 
is now depressed and this sends extra current thru the B coils (but not 
enough to pick up any of the relays). The operate switch is again set to 
position 3, and only the B coils now receive any current. (Any non-B 
relays that had previously been operated drop out, leaving only 
o(A, B,C, D) & a & B impressed on the machine.) At this point 
Mays states “the Ad pointers fall to the vertical.’’ I assume this is 
intended to mean, when o(A, B, C, D) isa tautology, that is, the initially 
cleared state of the machine, that only the Ab pointers fall to the ver- 
tical, that is to say, all the other relays remain energized. Under the 
present interpretation the AB and the Ab pointers (that is, the A 
pointers) fall on the first step and the ab pointers fall on the second 
step. This leaves only the aB pointers still horizontal. 

If the premise a V B is to be recorded, we must not drop cither the 
AB or ab pointers. Under the present interpretation this can be done 
in one step. With o(A, B, C, D) stored in the machine, we have the 
operate switch resting on 2. Literal keys A and 6 are simultaneously 
depressed. This sends additional current through the a coils and also 
through the B coils (but not enough to pick up any of the relays, even 
the aB relays, which get a double increment of current). The operate 
switch is then set to position 3, and the only coils receiving current now 
are the a coils and the B coils. All relays corresponding to atoms that 
imply av B (= ~ (Abd)) that had previously dropped out remain de- 
energized, and relays corresponding to atoms that do not imply av B 
drop out now leaving 


B, C, D) & (av B) 


impressed on the machine. The operate switch is then set to position 2 
and the above result is stored. 

The principal objection to Mays’s analysis is that there seems to 
be no distinction between the & and V operators. For instance, 
suppose we wish to impress ~ (Ab) & ~ (BcD) which is equivalent to 
(av B)& (bv Cvd). In the present analysis the A and } keys are 
depressed simultaneously, the operate key manipulated, the A and } 
keys released, then the B and c and D keys are depressed simultaneously ; 
in other words each premiss is handled as a unit. On the other hand, 
Mays says that each literal is recorded separately, even within a premiss. 
When it comes to starting a new premiss there is no clear specification 
of how to start. He says ‘“‘We must thus continue to impress as many 
of the premisses of a problem as we please. .”’ But he has just been 
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talking about individual literals of a premiss, and if we continue in the 
same way we are just augmenting the first premiss with more literals. 
Furthermore, since he manipulates the operation switch between literals, 
he is, under the present physical assumptions, considering each literal 
to be a separate premiss. This means that ~(Ab) instead of being 
entered as a V B, is entered as a & B, that is, as ~ A & ~ 5, and the 
De Morgan law is violated. 

In conclusion certain simplifications and generalizations of the 
machine may be described. 

1. The battery B, is dispensable if we substitute a mechanical device 
to close the relay armatures, each of which will remain closed with the 
operate switch on 2. Return all coils and contact 2 to battery B». 

2. Neither Marquand nor Mays mentions that, owing to duality, 
the machine can be used to express truth functions in (not necessarily 
distinguished) disjunctive as well as conjunctive normal form. Let 
operated relays correspond to falsity (redesignating the fiducial state), 
capital letters to true variables, simultaneous depression of literal keys 
to conjunction, and the use of the operate key to disjunction. The 
effect now is to accumulate true atoms, by losing falsities as the relays 
drop out. This is a more natural mode of operation, if we wish to com- 
pare the machine with present-day logic machines. 

3. The machine generali zesnaturally to n variables, with the limita- 
tion that the amount of equipment approximately doubles when n goes 
ton+ 1.7 Also, 7/t > n+ 1. 


7 This limitation is not inherent. See WiLtiAM Mieace, “Burroughs Truth-Function 
Evaluator,” J. Assn. Comp. Machinery, Vol. 4, pp. 189-192 (1957), and Artnur W. Burks, 
Don W. WarRREN, Jesse B. Wricut, “An Analysis of a Logical Machine Using Parenthesis- 
Free Notation,"” Math. Tables Aid. Comp., Vol. 8, pp. 53-57 (April 1954). 
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LIFETIME OF THE 161-KEV EXCITED STATE IN Sn’ * 


BY 


F. R. METZGER 


ABSTRACT 


Using the thermal resonance fluorescence method, the mean life of the 161-kev 
M1 transition in Sn"? has been measured as (5.77}3) X 10~" second. In arriving at 
this value, a total conversion coefficient 8r = 0.15 was assumed. 


INTRODUCTION 


The 161-kev transition in Sn''’ is known to be almost pure M1, 
the reported (1)! mixing ratio being E2/M1 = 0.001570:0013. It is, 
therefore, not surprising that Fagg et al. (2), using 4-Mev alpha par- 
ticles, were unable to observe Coulomb excitation of the 161-kev level 
in spite of the use of enriched Sn"'’. Since no information concerning 
the lifetime of this /-forbidden M1 transition was available, a measure- 
ment using the resonance fluorescence technique was undertaken. Both 
the thermal and the centrifuge methods are suitable for this isotope. 
For the available source strength and source purity, the thermal method 
was judged to be capable of giving a more accurate result for the mean life. 

The angular distribution of the resonance radiation for a spin se- 
quence 4 — 3 — 3 is of the form W(@) = 1 + AP, (cos@), that is, 
does not contain higher terms for any E2/M1 mixture. The total 
scattering cross section may, therefore, be obtained from a single 
measurement at a scattering angle (@ ~ 125°) for which P, vanishes. 


EXPERIMENTAL PROCEDURE 
A. General 


When the same isotope is used as the source of the exciting radiation 
and as the scattering material, the resonance condition is usually, at 
least partially, destroyed by the recoil energy losses in the emission 
and absorption processes. The thermal method (3) is one way of 
restoring the resonance condition. In analyzing the results obtained 
with the thermal method one assumes that the temperature dependence 
of the overlap of emission and absorption line is that calculated on the 
basis of Maxwellian velocity distributions in source and scatterer. A 
measure of this overlap is the average cross section (averaged over the 


__ * Assisted by the U.S. Office of Naval Research. Reproduction in whole or in part is per- 
mitted for any purpose of the United States Government. 
' The boldface numbers in parentheses refer to the references appended to this paper. 
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incident line) which may be written in the form (3) 


Cay = 2k(T; ra T:) X] 2k(7T, + 


where g, and g, are the statistical weights of excited state and ground 
state, respectively, \ is the wavelength of the gamma radiation of 
energy E,, T is the total width of the level, [') the partial gamma width 
of the transition to the ground state, M the mass of the recoiling nucleus, 
k Boltzmann’s constant, ¢ the velocity of light, 7, the absolute tem- 
perature of the source, T, the absolute temperature (4) of the scatterer. 

The presence of non-resonant elastic scattering, of room background 
and of contributions from other transitions makes it difficult to determine 
gay on the basis of a single measurement. It is much easier to measure 
the change in o,4, brought about by a change in the source temperature 
and to evaluate I'»?/T from this change of ¢4,. 


B. Source Preparation 


Ten milligrams of SnO:, enriched? to 92.62 per cent in Sn''*, were 
irradiated for a period of one week in the Oak Ridge National Labora- 
tory Research Reactor in a slow neutron flux of 1.7 & 10" n/cm? sec. 
In this way approximately 0.5 millicuries of 14-day Sn''? were produced. 
Impurities giving rise to gamma rays of higher energies were found to 
represent less than one per cent of the initial activity. They did, 
however, through Compton scattering, contribute appreciably to the 
non-resonant background counting rate in the pulse height region 
corresponding to the 161-kev radiation. For the studies reported here, 
no effort was made to remove these impurities, especially since an 
appreciable fraction of the hard gamma radiation originated from 
118-day Sn'", 


C. Geometry 


The experimental arrangement used for the scattering studies is 
shown in Fig. 1. The cylindrical tin scatterer was 8.3 in. in diameter, 
2 in. high, and 0.080 in. thick. Cd was used for the comparison scat- 
terer. Measurements were taken with the source at room temperature 
and with the source at 700°C. The counting rates were of the order of 
1000 counts per minute in a channel comprising approximately the full 
width at half maximum of the 161-kev line. The increase in the scatter- 
ing from the tin scatterer upon heating of the source from room tem- 
perature to 700°C. amounted to approximately one per cent of the total 
counting rate. The cross sections o,4,(700°C.) and o,,(27°C.) were 


2 Obtained from the Stable Isotopes Division of the AEC. 
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determined absolutely, that is, not by comparison with the Rayleigh 
cross section. The absolute determination was chosen because it was 
felt that the identification of the non-resonant counting rate in the 
elastic peak with the Rayleigh cross section is highly problematical even 
for a source emitting a single gamma ray. If, as in the case of the Sn'"” 
source used in this experiment, higher energy gamma rays are present, 
this identification would clearly be incorrect. In order to obtain the 


Oven 


_- Source 


Au 
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Scatterer 


0.016 in. Al + 0.008 in. Pb -+ 


Nal (TL) Crystal 


inches Photo - 


multiplier sin. Pb 
RCA 6342 


Fic. 1. Experimental arrangement. The temperature of the source was measured with a 
Chrome! P-Alumel thermocouple. 


absolute cross sections, the counting rate with the detector, a 4} in. 
thick, 1} in. diameter Nal (T1) crystal, exposed to the direct beam from 
the unshielded source was compared with the resonance counting rate, 
taking into account the geometry of Fig. 1. Since the exciting radia- 
tion and the resonantly scattered gamma rays are essentially of the 
same energy, this comparison does not involve the knowledge of the 
efficiency of the detector for different gamma ray energies. 


RESULTS AND DISCUSSION 


With 7, + 7, = 610° for the experiment carried out at room tem- 
perature, and 7, + 7, = 1281° for the high temperature experiment, 
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one calculates from Eq. 1: 


Agay = oay(700°C.) — o4,(27°C.) = 1.38 K 10-9 cm? (2) 


where I’, and I are expressed in ev. 
Fer the geometry depicted in Fig. 1, the resonance counting rate 
difference is calculated as 


ANres = 1.95 4 10'°N,AAcay 


(3) 


where Ny = C(d) X d? is the product of the direct beam counting rate 
and the square of the distance source-detector (in inches) employed in 
the direct beam measurement, and A is the natural abundance (0.0761) 
of Sn''?. Before calculating No from the direct beam counting rate 
C(d), the measured C(d) was reduced by 3 per cent. This reduction 
represents the estimated contribution of the gamma rays from the 
highly converted 159-kev M4 transition preceding the 161-kev M1 
radiation of interest. In arriving at this contribution, theoretical con- 
version coefficients (5) for pure M4 and pure M1 transitions were used. 
Combining Eqs. 2 and 3, one has 


= 4.88ANees/ No. 


(4) 


With the theoretical (5) total conversion coefficient for a pure M1 transi- 
tion, Br = 0.15, [,/T takes on the value T/T! = 0.87, and Eq. 4 goes 
over into 


= Es No. (5) 


From ten separate runs, an average value of (20.6 + 7.7) & 10-* was 
obtained for the ratio ANges/No. With Eq. 5 this then leads to a 
partial width IT) = (1.16 + 0.43) K 10-* ev, or a mean life 


tr, = (5.7733) 10-" second 


for the 161-kev transition to the ground state of Sn". With [,/T 
= 0.87, the mean life of the 161-kev level becomes ryye = (5.047%) 
second. 

The partial M1 lifetime is approximately two orders of magnitude 
longer than the single particle estimate (6). This is, however, in line 
with other M1 transitions and does not indicate an effect of the /-for- 
bidden character of the transition. The partial E2 mean life is longer 
than the single particle estimate or at best equal to it. 
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AUTOMATED COORDINATOGRAPH 


BY 


C. W. HARGENS 


Many uses will be found for the new automated coordinatograph 
developed for Aero Service Corporation, Philadelphia. The develop- 
ment project has transformed the respected Haag-Streit coordinatograph 
into a completely automatic plotting and measuring instrument. The 
aim achieved in this instance was to automate the instrument without 
losing the 0.001-in. accuracy obtainable in its original, manual form. 

Layouts on metal or glass with this degree of accuracy and nearly 
four feet in extent are now made from pre-programmed data. This is 
the plotting function. 

In the reverse sense tabulations of coordinate values may be rapidly 
produced from aerial photographs, astronomical plates, templates, and 
any other precise, two-dimensional representation. The dimensional 
values are at the same time punched on tape so that entry into a com- 
puter or other data processing arrangement is simplified. 

The equipment (Fig. 1) consists of the coordinatograph table, an 
electronic console, and the input-output unit. The control equipment 
was designed so that it could be added to the Swiss instrument without 
drilling holes or making similar alterations. This permits those who 
already have instruments to add the numerical control components 
without extensive modification. In other respects the quality of the 
electro-mechanical drives was made compatible with the original 
mechanism. ‘The scheme is remarkably foolproof in operation, as is 
indicated by the absence of limit switches, etc. 

The movement of the X- and Y-axis carriages is produced by means 
of a novel synthesis of steel bands and short leadscrews. This com- 
bination makes it possible to move rapidly over large distances and yet 
locate precisely to one thousandth of an inch. It is possible also to 
change driving speeds as much as 10,000 to one. A differential and 
electric clamp form a part of this mechanism. In the automated 
instrument speed changing is done automatically, a decade at a time, 
as the more significant digits reach coincidence with the numerical 
commands. 

The drive system of the new coordinatograph is completely divorced 
from the distance-measuring or feedback parts. The original rack and 
pinion which provided the accuracy of the Haag-Streit instrument are 
retained for their intended purpose. The reading dials have been re- 
placed with highly-resolved digital encoders. These units are capable 
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Fic. 1. Photograph of equipment. 


of reading five decimal digits and therefore can indicate the extremes 
of carriage motion in either the positive or negative direction. 

It will be appreciated that five decimal digits require at least 20 
binary digits. There are actually more than this because of lead-lag 
circuits associated with the most significant digits of the encoders. 
The conductor system of each moving carriage has nearly 60 circuits 
including the motor provisions. The digital data are transmitted from 
the position encoders by means of flat, Mylar sandwich conductors 
without any sliding contacts or slip rings. 

The block diagram (Fig. 2) shows the general arrangement of the 
feedback and control circuits. The binary-coded, cyclic decimal data 
from the coordinatograph are translated to Arabic decimal through 
transistor-buffered relay trees and odd-even complementing circuits. 
Translations and storage functions are obtained by relay arrays; 
numerical comparisons and decisions are made by means of diode 
matrices. There are no vacuum tubes in the equipment and no analog 
circuits to drift. Pulse counting is avoided on the grounds that a 
static logic system, directly coupled at all times, is more reliable. 

Operation of the comparator circuits is such that decimal digits are 
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2. 


examined sequentially. The decimal value of the command is available 
from the tape translator-storage unit while information in decimal form 
describing the carriages’ location is also sent to the comparator. The 
most significant digit is compared first, and a ‘‘decrease’’ or ‘‘increase”’ 
circuit is closed to the high-speed band motor. When coincidence be- 
tween the command and the carriage location is reached, the next most 
significant digit’s comparator circuit is enabled to function. This 
process continues, using various drive combinations, until complete 
coincidence of all X and Y digits is reached. At this instant all motors 
are stopped, and a prick mark is made. 

If the instrument is set for automatic plotting, the memory will 
be cleared and new data from the tape reader will be entered for the 
next cycle. The manual controls may be operated at any time without 
interference with the tape commands. Manual control circuits merely 
suspend memory control temporarily. When it is desired to again 
place the memory in control, this may be done without loss of the 
original entry. 

At any time the instrument will print out on demand the carriages’ 
X-—Y locations as shown on the numerical displays at the top of the 
console. If the operator wishes, he may also enter commands directly 
from the electric typewriter keyboard, without first preparing a tape. 

The equipment is designed on a building block principle so that 
special modifications can be made available for particular applications. 
A unit for “reading only”’ is such an example. More elaborate systems 
include line-following and closed-circuit television scanning equipment. 
A number of marketable units are now under construction by the sponsor, 
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BOOK REVIEWS 


THEORY OF RELaATiVviITy, by W. Pauli, trans- 
lated from the German by G. Field. 241 
pages, diagrams, 6 K 94 in. New York, 
Pergamon Press, 1959. Price, $6.00. 


Almost forty years have passed since Pauli 
first published (in German) this work on rela- 
tivity in ‘“‘Encyklopddie der Mathematischen 
Wissenshaften”’ (Oct. 19 Leipzig, 1921). An 
enormous literature on relativity has been 
published since. Yet Pauli’s treatment remains 
probably the most lucid and direct treatment 
of this subject. The evolution of ideas (from 
Lorentz to Poincaré to Einstein) that led to 
the postulates of special relativity is carefully 
described. The fascinating historical remarks 
throughout the book, showing the hesitations 
and misconceptions of the “founders” con- 
tribute to improve the understanding of this 
subject matter. About forty pages are devoted 
to general relativity. The last thirty pages 
describe various theories of elementary parti- 
cles (Weyl's, Mie’s, Einstein's) 

The translation of the original German text 
is followed by a substantial section of notes 
written by Pauli in 1956, especially for this 
edition. 

HENRI AMAR 
The Franklin Institute Laboratories 
and Temple University 


THERMODYNAMICS, by Herbert B. Callen. 376 
pages, diagrams, 6 X 9 in. New York, 
John Wiley & Sons, Inc., 1960. Price, $8.75. 


This text, written for advanced undergradu- 
ate and first-year graduate students, develops 
a logical theory of thermodynamics by logical 
deduction from a number of basic postulates. 
This abstract formulation is shown to result 
in a great simplification of the mathematics, 
and helps the reader appreciate the logical 
consistency and basic simplicity of the 
subject. 

In Part 1, comprising about half of the book, 
the general principles of classical thermody- 


namics are treated with the help of basic con- 


cepts and postulates, and some more advanced 


mathématical formalism which is supple- 
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mented by brief summaries in the appendices. 
Subjects include equilibrium conditions, proc- 
esses and thermodynamic engines, system 
stability, and phase transitions. 

Part II illustrates representative applica- 
tions of the method to chemical thermody- 
namics, electric behavior of solid systems, and 
magnetic and electric systems. 

Part II] demonstrates the advantages of 
the postulational approach in the study of 
fluctuations and irreversible thermodynamics, 
including thermoelectric and thermomagnetic 
effects. It is shown that the division between 
conventional thermodynamics and sstatisti- 
cal mechanics may be shifted so as to in- 
clude these subjects within the realm of 
thermodynamics. 

As the author states in his Preface, 
“. . . thermodynamics is the last of the major 
classical theories to undergo a postulatory rev- 
olution."’ The advances made possible in me- 
chanics by Lagrange’s and Hamilton's formu- 
lations, and in electromagnetic theory by 
Maxwell's equations, justify the hope that 
the endeavor towards more inclusive abstrac- 
tions, such as the one presented in Professor 
Callen'’s book, will be equally fruitful in the 
field of thermodynamics. The writing is ex- 
ceptionally clear, and retains a concrete flavor 
even in the presentation of the most abstract 
concepts and interpretation of mathematical 
relationships. Both as a text and a reference, 
this book may be highly recommended. 

The Franklin Institute Laboratories 


Vistas IN Astronomy, Volume 3, edited by 
A. Beer. 345 pages, illustrations, 7} K 9} 
in. New York, 1960. 
Price, $18.00. 


The third in what is called a permanent 
series of astronomical vistas is now off the 
press and, while it is a smaller and more 
specialized volume than its predecessors, it 
maintains the scientific impact of the initial 


Pergamon Press, 


volume. 
The international character of astronomy 
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in the book is once more demonstrated by 
having 21 authors representing six countries 
contribute articles to this volume. More than 
half the authors are British with the majority 
of the remainder working in the United States. 

Vistas in Astronomy, Volume 3 is divided 
into two sections. The first—in astronautical 
jargon—is on near space, that is, astronomy 
to the limit of the solar system. The second 
section encompasses the stars and galaxies. 

In the first section is found an extreme 
range of topics. These deal with the redis- 
covery of Neptune to the design of larger 
telescopes to the behavior of the recent comet- 
ary visitors. 

To this reviewer the sections on the Clock 
Paradox are fully and skillfully handled by 
Professors W. Cockran of Cambridge, Eng- 
land, and V. Fock of Leningrad, U.S.S.R. 
Academician Fock points out quite clearly 
what others neglect to do in the discussion 
of the Clock Paradox and that is that one 
body receives an acceleration and, thus, the 
two frames of reference are not equivalent. 

The last four contributions are on the two 
naked eye comets of recent years. Contained 
in the Porter contribution is the explanation 
for the cometary ‘‘beard"’ of the Arend-Ro- 
land comet. 

Part II of the book is a noteworthy con- 
tribution to astronomy. Prof. John B. Irwin 
presents a challenging picture as to what the 
astronomer could do were he given basic data 
on 70,000 stars. He develops the thesis that 
our basic knowledge of the brighter stars is 
indeed limited and this must be improved. 
He indicates that by using the photomultiplier 
tube in place of the photographic plate, it is 
possible to reach the 70,000 stars with a new 
order of accuracy. Prof. Irwin indicates that 
four observatories distributed in latitude in 
four to five years could complete this survey. 
Enough basic data could be accumulated to 
provide a breakthrough in our knowledge of 
stellar evolution. Prof. Irwin presents a sig- 
nificant argument and it is hoped that this 
will be implemented in the near future. 

Dr. Olin J. Eggen has contributed one of 
the most significant articles in the book. Like 
Dr. Irwin, he is concerned with stellar evolu- 
tion and he would like to attack the problem 
from the color-luminosity plots of the stars. 
Dr. Eggen has undertaken a most challenging 
problem and this article may provide a spring- 
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board for subsequent studies which may throw 
a good bit of light on the evolutionary history 
of the stars. 

Dr. L. Mestel discusses the role of the ga- 
lactic magnetic field in the condensation of 
interstellar matter which may eventually be- 


come nuclei for stars. 
The observational aspects of cosmology are 
explored by Dr. D. W. Sciama as the final 


contribution. 

The editors and publishers of this series are 
to be congratulated in their assembling of the 
vast talents represented between the covers 
of Vistas in Astronomy. 

I. M. Levitt 
The Fels Planetarium 


A BeGinner’s GUIDE TO THE Stars, by R. 
Newton Mayall and Margaret Mayall. 
184 pages, charts, 44 X 7} in. New York, 
G. P. Putnam’s Sons, 1960. Price, $2.50. 


Today astronomy is news. In tonight's sky 
you will find a complex of earth satellites and 
in the newspapers we find hints as to where 
in the sky to look to see them. With the 
advent of astronautics there has suddenly ap- 
peared a renewed and highly developed in- 
terest in the skies. Several books have ap- 
peared to aid those with this interest and 
guide them in learning the constellations. 

One such book is A Beginner's Guide to the 
Stars. Written by two highly qualified as- 
tronomers, the book is a magnificent primer 
to use in the introduction to the sky. 

The book is divided into two sections. The 
first deals with what is in the sky such as 
stars, moon, planets, comets, clusters and, 
finally, the man-made objects—the satellites. 
The second half of the book is an array of 
charts of the constellations with significant 
information concerning the star group. 

Accompanying the map are found the time 
of visibility of the constellation, its location, 
the characteristic identifying marks or fea- 
tures of the constellation and, finally, some 
notes in which we learn whether the constel- 
lation is a zodiacal, where are the interesting 
telescopic objects and other points of interest. 

The Mayalls have put together a most read- 
able pocket-sized book which is highly recom- 
mended for the beginner. 

1. M. Levitt 
The Fels Planetarium 
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ArROSPACE Dictionary, by Frank Gaynor. 260 pages, 6 X 9 in. New York, Philosophical 
Library, 1960. Price, $6.00 


DESIGN FOR A BRAIN, by W. Ross Ashby. Second edition, 286 pages, diagrams, 54 X 8} in. 
New York, John Wiley & Sons, Inc., 1960. Price, $6.50. 


EVOLUTION: Process AND Propuct, by Edward O. Dodson. 352 pages, illustrations, 6 X 9 
in. New York, Reinhold Publishing Corp., 1960. Price: $5.75 (college); $6.90 (trade). 


Fats AND OiLs, by H. G. Kirschenbauer. Second edition, 240 pages, diagrams, 6 X 9 in. 
New York, Reinhold Publishing Corp., 1960. Price, $7.00. 


SEMICONDUCTOR ABSTRACTS. VOLUME V—-1957 issue, edited by C.S. Peet. 449 pages, 84 X 11 
in. New York, John Wiley & Sons, Inc., 1960. Price, $12.00. 


FUNDAMENTALS OF SEMICONDUCTORS, by M. G. Scroggie. 160 pages, diagrams, 54 X 8} in. 
New York, Gernsback Library, Inc., 1960. Price, $2.95 (paper). 


Arps To Brocnemistry, by S. P. Datta and J. H. Ottaway. Fifth edition, 266 pages, diagrams, 
4 X 64 in. London, Balliere, Tindall and Cox, 1960. Price, $3.75 (available from The 
Williams & Wilkins Co., exclusive I 


. S. agents). 


THe Story or CHemistry, by Georg Lockemann. 277 pages, 44 X 7} in. New York, 
Philosophical Library, 1960. Price, $4.75. 


SCIENCE AND ENGINEERING IN AMERICAN INDUSTRY. Report on a 1956 Survey. 117 pages, 
8 X 10} in. Washington, U.S. Government Printing Office, 1959. Price, $0.70 (paper). 
NSF 59-50 


STATISTICAL HANDBOOK OF SCIENCE EDUCATION. 94 pages, 6 X 9 in. Washington, U. S. 
Government Printing Office, 1960. Price, $0.55. NSF 60-13 


Tue RestLess Atom, by Alfred Romer. 198 pages, diagrams, 4 7 in. New York, Double- 
day & Company, 1960. Price, $0.95 ($1.10 in Canada) (paper). Science Study 
Series 5 12 


Horns, Strincs & Harmony, by Arthur H. Benade. 271 pages, diagrams, 4 X 7 in. New 
York, Doubleday & Company, 1960. Price, $0.95 ($1.10 in Canada) (paper). Science 
Study Series 5 11. 


leuschler and R. Adler. 446 pages, diagrams, 6 X 9 in. 
Price, $12.00 


THE Soi AND Its Fertitity, by H 
New York, Reinhold Publishing Corp., 1960. 


BOOK NOTES 


SELECTED SemiconpucTOR Circurrs HANpD- information. A very wide range of subject 
ROOK, edited by Seymour Schwartz. V. p., matter is covered, such as: direct coupled 
diagrams, 6 X Yin. New York, John Wiley amplifiers, low-frequency amplifiers, high fre- 
& Sons, Inc., 1960. Price, $12.00 quency amplifiers; oscillators, switching cir- 


cuits, logic circuits, a-c to d-c power supplies’ 


Vhis practical handbook provides the elec- 


tronic engineer with a reference of over 150 converters ; small signal nonlinear circuits, and 


transistorized magnetic circuits. The latest 


usable, well designed, reliable semiconductor 
circuits representing the most recent advanced — type | R.E. standards are generally used to 


engineering techniques as well as design theory provide uniformity throughout the book. 


148 


| 
| 
2 
: 
: 
: 
: 
3 
: 


Aug., 1960.] Book 


METALLIZ1NG OF PLastics, by Harold Narcus. 
208 pages, illustrations, 5} K 74 in. New 
York, Reinhold Publishing Corp., 1960. 
Price, $5.50. 


Here is yet another addition to the Plastics 
Applications Series. Complete details are 
presented in this handy little book, for car- 
rying out every commercial metallizing pro- 
cess for plastics or other non-conductors. It 
is the first treatment of the subject to deal 
with actual production procedures, formula- 
tions and techniques for all known metallizing 
methods. 
cess (developed by the author), the deposi- 


These include the copper film pro- 


tion of “electroless” nickel coating, a review 
of the new molded conductive plastics, “gas"’ 
plating, the deposition of thick evaporated 
films and many others. Numerous illustra- 
tions are included showing the latest equip- 


More 
than a bibliographic source, this book is es- 


ment being used in the newer processes. 


sentially a metallizing manual for the plastics, 
electronics and electroplating industries. 


SoME MATHEMATICAL METHODS OF Puysics, 
by Gerald Goertzel and Nunzio Tralli. 300 
pages, diagrams, 54 & 8} in. New York, 
McGraw-Hill Book Co., Inc., 1960. 
$8.50. 


Price, 


This new work, designed primarily for phys- 
icists, presents a common set of mathematical 
techniques useful in such fields as quantum 
mechanics, acoustics, electromagnetic theory, 
reactor physics, etc. It is the outgrowth of 
material developed by the authors for gradu- 
They 
have sought to given on mathematicians a well- 


ate courses at New York University. 


rounded and thorough indication of the basic 
concepts involved in the study of linear op- 


erators, vector spaces, orthogonal functions 


and Fourier analysis, Green's functions as 


tools for the solution of physical problems. 


StyLe Guipe For Cuemists, by Louis F. 
Fieser and Mary Fieser. 116 pages, 6 X 9 
in. \ New York, Reinhold Publishing Corp., 


1960. Price, $2.95. 


This book considers the problems facing the 
writer of a thesis, report, or journal paper and 
presents helpful suggestions for recognition 
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and avoidance of errors. Positive advice is 
offered on writing techniques for enlivenment 
of style without aspiring to a literary level. 
A scientist can learn to put style into his 
sentences that will hold the reader's attention. 
Commentaries on points of grammar and rhet- 
oric, and suggestions on speaking are included 
as well as information on spelling, abbrevia- 
punctuation, choice of 


tions, hyphenation, 


words and Latin and Greek prefixes. 


PROBABILITY: AN INTRODUCTION, by Samuel 

Goldberg. 322 pages, diagrams, 8} X 5} 
New York, Prentice-Hall, Inc., 1960. 
$5.95 


in. 
Price: $7.95 (trade edition); 
edition). 


(text 


This new volume deals fully with the es- 
sentials needed by all who use probability con- 
cepts, apply statistical techniques, or measure 
the effects of chance phenomena. Logical in 
its approach, the book formulates its proba- 
bilistic theory by employing the language and 
notation of sets. By emphasizing experiments 
with a finite number of possible outcomes, the 
study requires no previous knowledge of cal- 
culus and can be readily understood by the 
reader with a good background in high school 
algebra 


Guass: Its INDUSTRIAL APPLICATIONS, by 
Charles J. Phillips. 252 pages, illustrations, 
6 X 9in. New York, Reinhold Publishing 
Corp., 1960. Price, $6.95. 


The manufacture, physical and chemical 
properties, and particularly, the applications 
of all types of glass, are covered in this vol- 
ume. It is directed toward engineers, scien- 
tists, production and purchasing personnel 
who are involved in the design, specification, 
testing, fabrication, purchase or sale of glass 
products. The author avoids both extreme de- 
tail and the popularized ‘‘no-data’’ approach 
by condensing a mass of information to tab- 
The book reflects all 


the recent advances in the field and includes 


ular or graphic form 


new information on such subjects as glass bead 
reflectorization and the “‘floot glass"’ process 
for making flat glass. Maximum continuity 
has been arranged so that each topic is pre- 


ceded by everything needed to understand it. 
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PHOTOCHEMISTRY IN THE LIQUID AND SOLID 
States, edited by Heidt, 
Robert S. Livingston, Eugene Rabinowitch 


Lawrence | 


and Farrington Daniels. 174 pages, dia- 
grams, 84 X 1llin. New York, John Wiley 
& Sons, Inc., 1960. Price, $6.00 


Based upon a symposium sponsored by the 
National Academy of Sciences and the Na- 
tional Research Council, this volume contains 
the contributions of outstanding authorities 
in the field of photochemistry. The authors 
present the basic principles of photochemistry 
storage, survey the field of photochemical re- 
actions, and state the requirements for re- 
action types which might prove useful for 
storing solar energy. Basic research findings 
are also presented with suggestions of areas 
for further research which can and will lead 
to the use of the sun as an important and 
inexpensive source of energy 


PrincipLes oF Unit Operations, by A. S. 
Foust, L. A. Wenzel, C. W. Clump, L. Maus 
and L. B. Andersen. 578 pages, diagrams, 
8} X11 in. New York, John Wiley & 
Sons, Inc., 1960. Price, $15.00. 


This book breaks away from the traditional 
treatment of unit operations concept and 
offers a new approach that is more economical 
in time. The unit-operations are presented 
as unified groups of operations stemming from 
identical fundamentals. The developments 
are built up from a simplified physical model 
or a basic mathematical relation, or both, 
using generalized notation. After a thorough 
coverage of the simplified models, the treat- 
ment progresses to the more complicated 
handling of realistic problems by applying the 


general equations to the specific operations 


AnaLytic Functions, by L. V. Ahlfors, H. 
Behnke, L. Bers, H. Grauert, M 
J. A. Jenkins, K. Kodaira, R 
and D. C. Spencer 
6 X 9 in. 
Press, 1960. 


Heins, 
Nevanlinna 
197 pages, diagrams, 
Princeton, Princeton University 
Price, $5.00 


This book contains one article each by the 
authors representing a comprehensive survey 


of recent developments both in the classical 


F. 1. 


and the modern fields of the theory of analytic 
functions. In the article by Ahlifors, the space 
of closed Riemann surfaces is shown to possess 
a complex analytic structure in which the 
periods of Abelian integrals are analytic func- 
tions. 
theory of complex spaces is presented by 
Behnke and Grauert. A proof of the Teich- 
muller-Ahlfors theorem based on the Beltrami 
Lindeloef’s prin- 


A survey of recent advances in the 


equation is given by Bers 
ciple and Martin's positive harmonic func- 
tions are the basic elements in a theory de- 
veloped by Heins. New and old results on 
the coefficient problems for Schlicht functions 
are derived by Jenkins by a unified method. 
Analytic functions on manifolds of two com- 
plex dimensions are treated in the article by 
Kodaira. 


proof of Liouville’s classical theorem on con- 


Nevanlinna presents a simplified 


formal mappings in Euclidean spaces of di- 
mension 2. Deformations of complex ana- 
lytic manifolds are studied by Spencer. 


THe Dynamic BEHAVIOR OF THERMOELEC- 
rric Devices, by Paul E. Gray. 136 pages, 
diagrams, 6 X 9in. Cambridge, The Tech- 
nology Press of M.1I.T.; New York, John 
Wiley & Sons, Inc.; 1960. Price, $3.50. 


This volume investigates the small-signal 
dynamic behavior of thermoelectric devices. 
Exact analysis of thermoelectric heat pumps 
and generators is difficult because these de- 
vices are described by differential equations 
that contain product-type nonlinearities. But 
since the dynamic behavior of these devices 
is of interest principally as a consequence of 
the requirement for control of the device, 
much useful information can be obtained by 
a small-signal analysis. This approach per- 
mits the treatment of linear models that des- 
cribe behavior of the devices. Analysis of 
these models produces a number of small- 
signal transfer-functions that can be used 
either to compute the response of devices in 
the frequency-domain, or to calculate the re- 
The 


results of this analysis can also be used to 


sponse of devices in the time domain. 


evaluate the quantitative effect of the ma- 
terials and device parameters on the dynamic 


behavior of devices. 
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INTRODUCTION TO STATISTICAL THERMODY- 
NAMics, by Terrell L. Hill. 508 pages, 
diagrams, 6 X 9 in. Reading (Mass), Ad- 
dison-Wesley Publishing Co., Inc., 1960. 
Price, $9.75. 


This book represents an introduction to 
equilibrium statistical mechanics, designed as 
a text for courses offered in chemistry or phys- 
ics departments on the advanced undergradu- 
ate-graduate level. The emphasis is on recent 
advances in statistical thermodynamics, and 
consequently the more traditional material is 
treated in a rather condensed manner to allow 
ample room for a survey of the important 
developments of the 1940's and 1950's. In 
keeping with the author's aim of providing 
only an introduction to the subject, topics are 
not treated exhaustively and at an advanced 
level. Thus, although the coverage of subjects 
is very broad, is it quite appropriate for a 
beginning text. 


CHEMICAL INSTRUMENTATION: A SYSTEMATIC 
APPROACH TO INSTRUMENTAL ANALYsIS, by 
H. A. Strobel. 653 pages, diagrams, 6 X 9 
in. Reading (Mass), Addison-Wesley Pub- 
lishing Co., Inc., 1960. Price, $9.75. 


Written at the advanced undergraduate- 
graduate level, this book may be used success- 
fully as a text for a course in instrumental 
analysis to two semesters in length. The 
book’s most distinctive feature is its attempt 
to convey a basic appreciation for the design 
and operation of instruments as a part of the 
larger scheme of making physical measure- 
ments. To facilitate this appreciation, back- 
ground chapters are presented on the science 
of instrumentation, error analysis, interaction 
of electromagnetic radiation and matter, phys- 
ical optics, electrical electrochemical theory 
Since the author's objective 
is to develop in the student a genuine under- 
standing of chemical instrumentation, he 
treats basic techniques fully. A good balance 
is maintained between physical theory and 
design on one hand, and chemical theory and 
procedure on the other. Above all, the text 
offers a systematic approach to the basic 
principles common to all methods and types 
of instruments, rather than a mere catalog 
of equipment which is rapidly outdated. 


and electronics. 
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CALCULUS OF FUNCTIONS OF ONE ARGUMENT, 
WITH ANALYTIC GEOMETRY AND DIFFER- 
ENTIAL Equations, by Edward J. Cogan, 
Robert Z. Norman and Gerald L. Thomp- 

587 pages, diagrams, 6 X 9 in. New 
Prentice-Hall, 1960. Price, 


son. 
York, 
$8.50. 


Inc., 


This study of calculus, completely in accord 
with modern views, applies the precise lan- 
guage of mathematics to the exploration of 
Without de- 
parting from traditional notation, this book 
breaks with classical treatment to promote 


concepts to their deeper roots. 


reasoning rather than mere problem solving, 
and emphasizes the function concept through- 
Employing contemporary terminology, 
the reader is aided in understanding the diffi- 


out. 


cult and advanced material by the use of such 
unifying ideas as solution sets, motivation by 
intuition and application, and integration of 
with and analytical 
geometry. Because of its independent view- 
point and modern approach, it will be espe- 
cially useful to readers interested in business 


calculus trigonometry 


administration as well as those specializing 
in science and engineering. 


CHEMISTRY OF NATURAL GUMS AND MUCcIL- 
AGES AND SOME RELATED POLYSACCHAR- 
ipes, by F. Smith and R. Montgomery. 
627 pages, diagrams, 6 X 9 in., New York, 
Reinhold Publishing Corp., 1959. 
$18.00. 


Price, 


The purpose of this book is two-fold: to 
evaluate critically the progress made in this 
fascinating field, and to encourage new ideas 
for the industrial application of gums. Every- 
thing that is known of the chemistry of these 
versatile substances is presented in this com- 
prehensive new monograph; this includes the 
The 
treatment of the subject is exhaustive, and 
begins with gum function and origin, The 
gums are classified botanically before the text 


rapid strides made in the last few years. 


takes up their chemical and physical prop- 


erties. Because of the growing commercial 
importance and possibilities of these raw ma- 


this 


with the sound basis of information it needs 


terials, monograph provides industry 
for the successful development of gums and 


mucilages. 
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NATIONAL BUREAU OF STANDARDS NEWS 


IMPROVED METHOD FOR MEASURING ABRASION RESISTANCE 
OF COATINGS 


The National Bureau of Standards has made a number of improve- 
ments in its Abrasive Jet Method for measuring the abrasion resistance 
of protective coatings. These improvements' have resulted from con- 
tinued studies by A. G. Roberts, who originally developed the method 
in 1955. The changes have not altered the basic manner of operation 
but have significantly increased the speed, ease, and precision of the 
measurements. 

In the basic method,’? pressurized gas propels an abrasive powder 
from a vibrating storage chamber through a nozzle onto the test speci- 
men. Abrasive flow, pressure, distance, and angle are all closely con- 
trolled. Abrasion resistance is measured by the time required to 
abrade through the coating and is expressed in terms of a unit thickness. 
The first show of bare substrate is the end point and is clearly detected 
by an abrupt change in color or gloss. The high speed of the particles 
in the abrasive jet provides a rapid rate of abrasion that permits most 
coatings to be tested in a matter of seconds. The abrading conditions 
can be varied over a wide range to permit the testing of virtually any 
type of coating. 

The small scale of operation obtained with the equipment allows 
the economical use of continually fresh abrasive as well as a large number 
of replicate tests on a single specimen. Thus the use of a reference 
panel for calibration and standardization purposes is feasible and has 
become an important part of the improved procedure. The values 
recorded for the reference panel also provide a reliable means for reduc- 
ing data obtained under different conditions to a common basis for 
comparison. The coefficient of variation for measurements on the 
reference panel is less than two per cent. 

Close control of the amplitude of vibration is necessary to assure a 
uniform flow of abrasive and a constant rate of abrasion. Recent stud- 
ies have shown that the amplitude of vibration changes with the gradual 
drop in the level of abrasive as operation continues, because of a shift 


‘For further details, see “Improved NBS Abrasive Jet Method for Measuring Abrasion 
Resistance of Coatings,” by A. G. Roberts, ASTM Bul., No. 244, p. 48 (1960). 

? “Abrasive Jet Method for Measuring Resistance of Organic Coatings,”” by A. G. Roberts, 
W. A. Crouse, and R. S. Pizer, ASTM Bul., No. 208, p. 36 (1955); and NBS Tech. News Bul., 
Vol. 39, p. 33 (1955). Apparatus for Measuring Abrasion Resistance, U. S. Patent No. 
2,907,200, Oct. 1959. 
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in the center of gravity of the electromechanical vibrator system. 
Since the change is gradual, it can be compensated by periodic readjust- 
ment of the regulated line voltage applied to the vibrator, in order to 
maintain the standard abrasion rate on the reference panel. Unde- 
sirable turbulence and interference effects in the abrasive storage 
chamber have been minimized by redesigning the chamber so that the 
abrasive particles sift into and are propelled along a straight tube 
instead of through a circular chamber. Use of a more free-flowing 
abrasive (aluminum oxide) than was formerly employed also has im- 
proved the uniformity of operation. 

New devices for setting and adjusting nozzle-to-coating distance 
have been added with which distance is quickly and easily set to within 
0.002 in. Several subtle sources of error in the original method have 
been identified and eliminated. Practical operating limits have been 
established within which good precision can be expected over a wide 
range of test conditions. 

Studies have shown that the method is not limited to the organic 
coating materials for which it was originally designed. It is also ap- 
plicable to ceramic, porcelain, and anodic coatings, as well as films and 
tapes when attached to a rigid substrate. Sheet metal of aluminum, 
brass, or steel may also be tested ; their rates of abrasion are of the same 
order of magnitude as those of the more resistant types of organic 
coatings. A good qualitative correlation has been observed between 
test results and the service performance of protective coatings applied 
to military aircraft. 

The speed and simplicity with which the measurements are made, 
the reproducibility of results, and the opportunity for calibration and 
standardization afforded by employing a practical reference panel are 
features of the improved method that offer a promise of considerable 
usefulness in the field of abrasion testing. 


SPECTRORADIOMETER FOR FLUORESCENT LAMPS 


A number of methods are available for evaluating the color of 
fluorescent lamps, and for studying the effects of the many variables in 
their manufacture. However, the expanding applications of fluorescent 
illumination today show that when these lamps are supplied by differ- 
ent manufacturers it is not always possible to achieve a perfect color 
match. The National Bureau of Standards has therefore constructed 
a spectroradiometer to calibrate fluorescent lamp color standards that 
may be used throughout the industry to obtain color uniformity.’ 
Measurements made with this instrument will also form a basis for 
international comparisons of fluorescent lamp colors. 


*For further technical details, see ‘“Ratio-Recording Spectroradiometer,” by H. K. 
Hammond, W. L. Holford, and M. L. Kuder, J. Research NBS, Vol. 64C, p. 151 (1960). 
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Spectroradiometers generally have three basic components: an 
clement for dispersing radiant flux from a light source into a spectrum, 
a radiometric detector for measuring this flux in each wavelength, and 
a recording device for plotting the resulting spectral irradiance data. 
From these data the color of the source is determined. Some spectro- 
radiometers‘ measure the current from the detector to find the spectral 
irradiance of fluorescent lamps. In the present work, however, a 
method was desired that would be essentially independent of both the 
response characteristics of the detector and the spectral transmittance 
of the dispersing element. A system was therefore designed by H. K. 
Hammond, W. L. Holford, and M. L. Kuder of the Bureau staff, in 
which a test lamp is compared directly with an incandescent lamp of 
known color temperature. 

A double-prism monochromator with a cam-linearized wavelength 
drive from 360 to 760 my is used as the dispersing element in the NBS 
instrument. The detector is a multiplier phototube with high spectral 
sensitivity throughout the visible spectrum. A 10-in. strip-chart re- 
corder is used to plot the relative spectral irradiance ratio of the two 
light sources. These three components, together with the power 
supplies and electronic circuitry, are housed in a console. The two 
light sources are mounted on the outside of the console at adjacent 
entrance ports. 

Optical attenuators, in the form of adjustable apertures, are located 
in the path of each light beam. Each attenuator has two identical 
components, one fixed and one movable. Each component is a flat 
metal foil with an array of 25 diamond-shaped holes, 1/10 in. ona side. 
The two components of an attenuator are fitted together in such a way 
that they slide over each other to enlarge or decrease the size of the 
holes as required. 

In operation, the spectroradiometer compares the radiant flux from 
each source as the spectrum is scanned. When the intensities of the two 
beams are unequal, the unequal phototube signals cause a servo motor 
to adjust the attenuators, simultaneously opening one and closing the 
other, until the beam intensities are equal. The attenuation required 
at each wavelength for equalization is automatically plotted as a func- 
tion of wavelength over the visible spectrum. The curve thus recorded 
is the ratio of relative spectral irradiance from the two sources. 

Wavelength correlation between monochromator and recorder is 
maintained by driving the prisms and the strip chart with a second 
motor. The speed of this motor is varied according to the magnitude 
of the unbalance signal. When the curve is flat, the rate of scan is 
maximum; when the curve becomes steep, the rate of scan is auto- 

*“An Automatic Recording Spectroradiometer,”” by V. K. Zworykin, J. Opt. Soc. Am., 
Vol. 29, p. 84 (1939); and “Spectroradiometry,”’ by F. J. Studer and W. R. Jacobsen, GE 
Review, Vol. 52, No. 10 (1949), 
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matically reduced to assure faithful recording of the peaks of sharp 
spectral lines. Auxiliary gears and clutches provide wavelength scales 
of either 1 mu or 10 my for each inch of chart length. At the latter 
rate, the time required for scanning the visible spectrum is about 8 
minutes. 

Although the spectroradiometer was built primarily for the calibra- 
tion of fluorescent lamps, it can also measure the relative spectral ir- 
radiance of any source or source-filter combination. Luminescent or 
fluorescent materials can be measured if they have several square inches 
of uniform luminance and if the luminance is at least one footlambert. 


YOU CAN ADVANCE SCIENCE EDUCATION 


Today, more than ever before in its 135-year history, there is vital need for 
The Franklin Institute to effectively promote education in science and technology. 
It is imperative that we meet this challenge by providing adequate educational op- 
portunities in these fields. This requires vision, objective planning, and money. 
We have more than enough of the first two requisites, but far too little of the third. 

Our programs are aimed at professional scientists and industry, as well as the lay 
public and young people seeking inspiration and guidance in choosing a career. 
The Institute’s educational programs are impressive, for they begin with students in 
the early grades of our elementary schools and continue throughout an individual's 
professional or industrial life. With more funds at our disposal, the scope and 
vigor of these activities could be greatly increased. 

The Franklin Institute is not richly endowed. It is a non-profit organization, 
depending for encouragement and support on an understanding public. Capable 
and conservative management assures wise administration of all funds. 

Your gift or bequest, large or small, will be deeply appreciated and will be used 
effectively to broaden the Institute’s educational usefulness. There is a warm 
satisfaction in giving financial support to an organization that has pioneered in, and 
is dedicated to, the advancement of science and technology through education. 

When property is transferred, title should be in the name of The Franklin 
Institute of the State of Pennsylvania for the Promotion of the Mechanic Arts. 

The Secretary of The Franklin Institute will gladly furnish you with additional 
information. Write to him at The Franklin Institute, Benjamin Franklin Parkway 
at Twentieth Street, Philadelphia 3, Pennsylvania. 
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Language Translation by Machine. 

The United States Air Force may 
have the answer to the age-old prob- 
lem of international tongue-twisting 
through its program of language trans- 
lation by machine. An electronic 
translator capable of turning Russian 
into English at the rate of 35 words 
a second has been developed under a 
five-year project at the Air Research 
and Development Command's Rome 
Air Development Center, Griffiss 


AFB, N. Y. 

The only one of its kind in existence, 
the translator has been in use on an 
experimental basis since June 1959. 
According to Air Force engineers, the 
English output is ‘‘no literary master- 
piece” since translation at this point 


is on a word-for-word basis. How- 
ever, the machine gives a very good 
idea of the general content of a Rus- 
sian article—even though it ignores 
many of the grammar rules. By the 
end of this year a sophisticated word 
analyzer with logic circuits for sen- 
tence structure will be added into the 
translation complex. 

Heart of the automatic language 
translator is a “‘photoscopic memory” 
invented by Dr. Gilbert W. King of 
International Business Machines Re- 
search Center, Yorktown Heights, N. 
Y. A transparent disk only 10 in. 
in diameter, this dictionary unit can 
store 550,000 Russian-English words 
in an area the size of a postcard. They 
appear around the edge of the disk in 
concentric tracks of binary code. As 
a Russian word is fed into the trans- 
lator via a punched tape, it is ‘‘read”’ 
by the machine and converted into 
electrical signals. These are matched 
in a lightning dictionary drill with the 
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coded equivalents on the glass disk. 
The English translation is transmitted 
over an electric typewriter. Any word 
in the dictionary disk can be located 
by the machine in less than 1/300 sec. 

The Intelligence Laboratory at the 
Rome Air Development Center is using 
the research and knowledge of top 
language scientists in this country and 
abroad. It has contracts with seven 
universities and four industrial firms 
in the automatic language translation 
program. 

The problem of machine translation 
falls into two essential parts: First, 
correlation of logic when basic lan- 
guage structures differ as much as 
they do in Russian and English; and, 
second, development of the machinery 
to relay this correlation into electron- 
ics. At present there are still some 
major obstacles in the way of fully- 
automatic, high quality translation. 

For instance, although translation 
can be completed at extremely high 
speeds—the translator scans 16,000 
words a second—the process of feeding 
material to the machine is compara- 
tively slow. 

The tape containing the Russian 
text is prepared manually by a typist 
working with a Russian keyboard 
typewriter. Using this method, only 
40 words a minute are prepared for 
translation. To overcome this lag, a 
special input machine is being designed 
for the Air Force by Baird Atomic 
Inc., Cambridge, Mass. A Russian 
print-reader, which will scan the ma- 
terial electrically and automatically 
make the required tape, will be ready 
for test later this year. It will bring 
the input capability up to 2400 words 
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a minute or sixty times faster than 
the human operator. 

A forerunner of the Russian print- 
reader which reads English electron- 
ically was successfully developed un- 
der contract to RADC by Intelligent 
Machines Research Company in May 
1959 and is presently in use. 

Another problem—that of transmit- 
ting graphic material along with the 
written word—is now under study for 
the Air Force by the Syracuse Uni- 
versity Department of Graphic Engi- 
neering. Design for the construction 
of a machine to handle this pictorial 
output will be initiated in May. 

The “thinking”’ portion of the trans- 
lator—the dictionary and word an- 
alyzer—is also being improved at IBM 
Research. A lexical ‘“‘buffer” has re- 
cently been added which stores words 
coming from the photoscopic diction- 
ary in complete sentences. The “‘an- 
alyzer” will provide corrected 
grammar. 

While the Air Force’s primary in- 
terest is in the translation of Russian 
technical material, the machine’s pres- 
ent vocabulary enables it to translate 
all Russian texts with the possible ex- 
ception of poetic language where the 
author may attach different meanings 
to certain words. If the translator 
does come across a word not on its 
glass listing, it transcribes it in red 
for later inclusion. 

The benefits to be gained for both 
scientific and academic purposes 
through the use of such a translator 
are many. At present the U. S. gov- 
ernment is spending over a million and 
a half dollars annually for the job of 
converting Russian publications into 
English. Because of the high expendi- 
ture of time and money required by 
human translation, only a part of the 
necessary work is being done. 

It is estimated that the Russians 
publish a billion words a year—one 
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third of these in scientific fields. The 
U.S. government translates 80 million 
words a year at an average cost of 
two cents a word. Even in its present 
state of development, the USAF auto- 
matic translator— which can be oper- 
ated by an individual with no knowl- 
edge of Russian—is able to render a 
page of Russian into English in half 
a minute compared to 45 minutes 
required by the average human 
translator. 

The Air Force language translation 
program, started in 1955, will be com- 
pleted late in 1961 when the experi- 
mental model of the fully integrated 
translation complex is unveiled. Esti- 
mates are that by that time five mil- 
lion dollars will have been spent on 
the project. Air Force engineers and 
scientists believe this amount is only 
a drop in the budget bucket compared 
to the returns that automatic language 
translation will pay off in the years 
to come. 


Turret System for Helicopters.-—.\ 
new lightweight turret system for 
helicopters which will permit accu- 
rate cockpit control of an externally 
mounted machine gun or surveillance 
TV camera has been developed by 
General Electric's Missile Production 
Section, Burlington, Vt., under a U. S. 
Army Ordance contract. The remote- 
control system is described as being 
easily and inexpensively adaptable to 
existing helicopters and light obser- 
vation aircraft. 

Capable of handling 150-lb. loads, 
it will accommodate any small machine 
gun up through 50 caliber to give 
low-speed aircraft protection against 
ground fire. It can also be used to 
provide remote-control for a conven- 
tional closed-circuit television camera 
or General Electric’s new Image Orthi- 
con “‘see-in-the-dark’’ TV camera for 
nighttime reconnaissance. 
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The system can be provided as a 
kit for installation by second-echelon 
field maintenance personnel, accord- 
ing to company engineers. Once in- 
itial installation is completed, it can 
be removed entirely or adapted in the 
field to either armament of surveil- 
lance equipment in a matter of min- 
utes to meet specific mission require- 
ments, they said. 

Previewed here at the annual meet- 
ing of the American Helicopter So- 
ciety, the system consists of a mount 
trainable in azimuth and elevation, 
drive motors, servo amplifiers, and a 
pistol-type operator hand control. It 
has no hydraulic components and uses 
standard 24-volt electric power. Total 
weight of the transistorized system is 
less than 80 pounds. 

The mount is capable of 370 degrees 
horizontal rotation and could provide 
a vertical range of plus or minus 45 
degrees. Quickly maneuverable, it 
provides slew speeds up to 70 degrees 
per second, the engineers said. 

The cockpit control, handled and 
aimed like a pistol, is suspended from 
the upper cabin framework by a paral- 
lelogram suspension system. ‘This sus- 
pension permits the co-pilot to freely 
move the pistol grip for clear aim from 
a natural position. The mechanism 
is stowed out of the way in the cabin 
when not in use. 

In addition to its airborne applica- 
tions, the system can be employed 
where unattended control of cameras 
or similar equipment in remote, inac- 
cessible or dangerous locations is de- 
sired. The engineers pointed to pos- 
sible future uses with closed-circuit 
TV cameras in plant or prison security 
systems, at missile launching sites and 
in radioactive atmospheres. 


Rough Terrain Crane.—.\ ‘‘sophisti- 
cated” 20-ton capacity rough terrain 
crane, designed for a variety of appli- 
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cations including operation in four feet 

of sea water, is now under test at the 
U. S. Army Engineer Research and 
Development Laboratories, Fort Bel- 
voir, Va. 

It is designed for use in military con- 
struction projects including clamshell 
work, dragline work, steel erection, 
bridge assembly, and other uses which 
are normal to crane operations. It 
also will be utilized for barge or landing 
operations, working out in the surf to 
assist in unloading landing craft. 

Weight of the basic 20-ton rough 
terrain crane, with a 30-ft. strut boom, 
is held to approximately 50,000 Ib. 
The crane is expected to be able to 
travel the highways, under its own 
power, at speeds of 30 mph., climb 40 
per cent grades at a speed of 2.25 
mph., negotiate 10 per cent side slopes 
and, by using maximum front and rear 
steering, turn inside a circle with a 
30-ft. radius. 

One of its unique features is a special 
system which permits the inflation and 
deflation of the four tires while the 
vehicle is standing still or in motion. 
Air pressure can be varied between 
10-85 lb., giving added traction when 
needed in certain types of soil con- 
ditions. Also, hydraulic motor-driven 
outriggers are permanently located on 
the carrier, and are lowered to help 
stabilize the crane when lifting maxi- 
mum loads. 

Mounted on a specially designed 
diesel-driven, off-highway type carrier, 
the crane is capable of operating with- 
out support equipment for sustained 
periods of time on and off the highway. 
It is even equipped with a dozer blade, 
with limited capabilities, to permit the 
machine to function independently of 
other construction equipment on the 
clearing and leveling of its own work 
site. 

The crane on the carrier is mounted 
on hook rollers and powered by a diesel 
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engine. Equipped with a 30-ft. strut 
boom or a 30-ft. gooseneck boom, and 
three 20-ft. insert sections, it is full 
revolving. Additional counterweights 
are provided for operation with a 90- 
ft. boom. Using the 30-ft. boom and 
outriggers, the unit can lift 10 tons 
at a 20-ft. operating radius. It has a 
25,000-Ib. drawbar pull. 

The prototype crane, currently un- 
dergoing tests, was built by LeTour- 
neau-Westinghouse Company, Peoria, 
IIl., under contract with the Labora- 
tories. The superstructure was man- 
ufactured by American Hoist and Der- 
rick Company. 


Improved Solid Lubricant.— <A solid 
lubricant with as much as 20 times 
the normal load-bearing capacity has 
been announced by Dr. Guy Suits, 
General Electric Vice President and 
Director of Research. Solid lubricants 
make possible the miniaturization of 
many devices, since they can be used 
at higher temperatures than the liquid 
type, Suits noted. They also provide 
permanent lubrication, improve relia- 
bility, and reduce maintenance. 

The improved lubricant was achieved 
by the addition of various inorganic 
sulfides to the well-known lubricant 
molybdenum disulfide, or to tungsten 
disulfide, Suits revealed. The effect 
of the additives was discovered by Dr. 
Arthur J. Haltner and Curtis S. Oliver 
of the General Electric Research Lab- 
oratory, and was evaluated in a long 
series of tests carried out in collabora- 
tion with other staff members. 

When a }-in. diameter chrome alloy 
steel hemisphere slid over a mild steel 
surface, using molybdenum disulfide 
as the lubricant, the lubricant film was 
able to withstand a one-quarter kilo- 
gram load. When an additive (silver 
sulfide) was incorporated in the lubri- 
cant, the load-bearing capacity in- 
creased to 5.3 kilograms. Even better 
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results are obtained when the lubri- 
cants are used on relatively hard sur- 
faces, such as chromium. 

The additives are all inorganic sul- 
fides, and their effects can be noted 
at relatively low concentrations. Dur- 
ing the tests conducted thus far, they 
were typically added in quantities of 
approximately 10 per cent by weight. 
Among the most effective of those that 
have been investigated are the sulfides 
of antimony, platinum, mercury, sil- 
ver, titanium, and lead. When used to 
lubricate iron or steel surfaces, such ad- 
ditives apparently decompose to pro- 
duce a layer of iron sulfide on the metal 
surface. This layer is believed to have 
better lubricating properties than the 
original metal, and may interact ad- 
vantageously with thé lubricant itself. 
The mechanism by which the addi- 
tives improve lubrication may be simi- 
lar to that which is effective when 
“extreme-pressure”’ additives are used 
in liquid lubricants. 

Studies carried out by General Elec- 
tric indicate that chemical reactions 
between molybdenum disulfide and the 
surrounding atmosphere may have a 
significant effect on lubricating effi- 
ciency, to an extent not previously ap- 
preciated. The investigations of mo- 
lybdenum disulfide are part of an 
extensive program being conducted at 
the GE Research Laboratory on the 
fundamental mechanisms involved in 
lubrication by solids. 


Inexpensive Method for Repro- 


ducing X-Ray Films.—A recent issue 
of The University of. Michigan Medical 
Bulletin reports that doctors and a pho- 
tographic specialist at the U-M have 
found a simple inexpensive method 
for reproducing X-ray films with “‘suf- 
ficient detail to maintain diagnostic 
quality.” Previous methods have 
been expensive and time-consuming, 
and have rarely produced satisfactory 
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copies. The Bulletin says the new de- 
velopment is highly important in get- 
ting maximum distribution of unusual 
X-ray pictures for medical consulta- 
tions and teaching purposes. 

The account was written by Dr. 
Walter M. Whitehouse (M.D.), asso- 
ciate professor of radiology ; Fred And- 
eregg, supervisor of U-M Photographic 
Services; and Dr. Fred J. Hodges 
(M.D.), chairman of the U-M Depart- 
ment of Radiology. 

They employed an improved diazo 
method of reproduction with a stan- 
dard diazo machine—somewhat like 
a blue-printer—used for reproducigng 
documents. 

However, the resulting duplicates 
differ from blue-prints in that they 
have a transparent base and a “‘direct 
positive continuous tone.”” A new 
type of diazo film produced by the 
Ozalid company made possible the ap- 
plication to X-rays. 

The process produces a dry finished 
copy in two minutes. It differs from 
the original only in having a thinner 
film base and a blue tint. 


Computer-Directed Navigation Sys- 
tem.—A computer-directed map pro- 
jection system under development by 
IBM engineers may enable the pilot 
of a supersonic plane to tell his loca- 
tion over the earth by merely glancing 
at a screen on his instrument panel. 

The experimental system can pro- 
ject a highly accurate circular map 
of an area 400 miles in diameter 
(125,000 square miles) on a 7}#-in. 
screen. An airborne computer per- 
forms all the calculations necessary to 
position the map on the screen and 
pinpoint the pilot's location. He can 
literally see on the display map what 
city or town he is flying over. 

Heart of the navigation display is 
a glass hemisphere about 6 in. in di- 
ameter. A detailed map of half the 
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Earth is reproduced photographically 
on its inside surface. A beam of light 
illuminates a small section of the 
map, which is then projected onto a 
flat, translucent screen in front of 
the pilot. As the plane moves, the 
computer automatically adjusts the 
map presentation. 

Engineers at the IBM Federal Sys- 
tems Division laboratory at Owego, 
N. Y., say that a similar display de- 
vice could be used in space navigation. 

To use the system in space, a tech- 
nique known as stellar map matching 
might be employed. A stellar map, 
or ‘chart of the heavens,” would re- 
place the map of the Earth inside the 
glass hemisphere and be projected on 
the screen. Then a telescope view 
from the space craft would also be 
projected on the screen, superimposed 
on the space map projection. When 
the stellar map and the telescope view 
are lined up, the computer can auto- 
matically determine the position of 
the space craft. 

IBM engineers are now working 
towards refinement of the system. 


Epitaxial Diffused Transistors. 
Major reductions in switching times 
and collector resistances of diffused 
base transistors have been realized at 
Bell Telephone Laboratories by the 
use of epitaxially grown films. These 
improvements have been achieved by 
combining diffused base technology 
with the epitaxial film technique. Films 
themselves have also been under in- 
vestigation in a number of other 
laboratories. 

Application of this technique has 
resulted in a reduction of the switch- 
ing time of silicon devices by a factor 
of more than 10, and in a comparable 
reduction in the collector resistance. 
This major development is expected 
to have far-reaching implications in 
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both the fabrication and application 
of semiconductor devices. 

Diffused base transistors require a 
relatively high resistivity collector re- 
gion in order to attain low capacitance 
and high voltage breakdown. Typi- 
cally, to afford ease of mechanical 
handling, this region has been much 
thicker than required electrically. The 
excess thickness increases the collector 
resistance and, through carrier stor- 
age, the switching time. 

Ideally the thickness of the lightly 
doped collector region should be in the 
neighborhood of 0.1 mil, which is about 
a factor of 30 thinner than normally 
used. Semiconductor wafers prepared 
by conventional methods become ex- 
tremely difficult, if not impossible, to 
handle as they approach this desired 
thinness. The problem can be solved 
by the use of epitaxial film techniques. 
Lightly doped epitaxial films are grown 
on and supported by a low resistivity 
substrate giving the desired combina- 
tion of electrical properties and me- 
chanical strength. The term “‘epitax- 
ial” means that the film is a direct 
extension of the single crystal struc- 
ture of the substrate wafer. 

Diffused base transistors have been 
made at Bell Telephone Laboratories 
on epitaxial layers of both germanium 
and silicon and have exhibited the 
theoretically predicted improvements. 
For example, in two similar silicon 
transistor structures, one conventional 
and the other using the epitaxial ma- 
terial, switching time in a typical cir- 
cuit has been reduced from 200 to 20 
millimicroseconds. Further, collector 
series resistance of the epitaxial tran- 
sistors was reduced by a factor of more 
than 10 and was comparable to that 
of conventional devices 15 times larger. 

Experiments with germanium indi- 
cate that use of epitaxial layers will 
extend the frequency response of ger- 
manium transistors well beyond that 
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of the two kilomegacycle device re- 
cently described by Bell Telephone 
Laboratories. 

The new technique involves modi- 
fication of the starting material and 
can be introduced directly into the 
front end of existing diffused base pro- 
duction facilities without any further 
changes. Single crystal wafers of 
heavily doped material are first cut 
and polished. A lightly doped thin 
film (about 0.1 mil) of the same con- 
ductivity type is then deposited epi- 
taxially on the wafer surface. This 
film provides the desired thin, lightly 
doped collector region. From this 
point on, standard diffused base tech- 
niques are employed. 

The use of epitaxial material in dif- 
fused base transistors not only results 
in major improvements in switching 
time and collector resistance, but in 
addition simplifies the design and un- 
derstanding of transistor devices and 
brings them closer to ideal forms, such 
as fi-p-i-n structures. Further, the ad- 
dition of the epitaxial film technique 
to the well established diffusion tech- 
nology provides the device engineer 
with an extra degree of design freedom 
which should result in new devices 
difficult or impossible to achieve by 
older techniques. 


Water Supply in Space.—-Research 
for a solution of the water supply 
for man on an extended space trip is 
underway at General Electric's Space 
Sciences Laboratory. The work is 
being done under a $94,000 contract 
recently awarded by the National 
Aeronautics and Space Administra- 
tion to the Company's Missile and 
Space Vehicle Department. 

Using a unique vacuum pyrolytic 
catalysis technique, the approach in- 
volves vaporizing the volatiles in the 
biological wastes at reduced pressures, 
hence low temperatures. The vapor 
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is then pyrolyzed over a catalyst and 
thence condensed back to a liquid 
state. This condensate has been an- 
alyzed and found to be without odor, 
color or toxic materials and adjudged 
suitable for human consumption, ac- 
cording to General Electric scientists. 
The technique does not require chem- 
ical additives or filtering. 

Rats and a monkey will be used in 
the study to determine of the feeding 
of recycled recovered water has any 
effect on animals. The rats will be 
studied through more than one gener- 
ation, while the monkey will be care- 
fully observed during a prolonged con- 
finement in a closed area. 

Besides visual observation, periodic 
tests of waste samples will be made 
to determine the animal's health. If 
indicated, the monkey will be re- 
moved for more thorough physical 
examination. 


CURRENT 


Topics F. 1. 

The monkey will be provided with 
a three-day supply of water in the 
semi-closed system. Dry food and 
dry air will be added as required. (In 
a closed system nothing is added or 
subtracted.) Waste products, includ- 
ing moisture from the vented air will 
be collected, processed, and returned 
to the system. 

The new approach on water recovery 
is being investigated owing to the 
weight and storage problems involved 
in supplying water by conventional 
methods, according to General Elec- 
tric scientists. 

The water supply problem is only 
one of many aspects of man’s environ- 
ment in space now being studied at 
the General Electric Space Sciences 
Laboratory. These investigations in- 


clude research on waste management, 
oxygen recovery, air purification and 
food requirements. 
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At a time when the focus is on outer space — a time when 
we know more about the moon's surface than we know 
about 70% of the earth’s surface — geophysicists such as 
this sea-going scientist from Texas Instruments are fur- 
— space-age technology with studies of inner space 
— the oceans. 


He represents one of more than 60 global TI geophysical 
investigation parties who are applying 30 years’ expe- 
rience in earth sciences to measurement of earth's land 
and sea physical characteristics. Fundamental to modern 
defense, earth sciences have extended their traditional role 
in petroleum exploration to gathering information vital 
to programs in underwater warfare, missile qecuracy, 
7 icatt , detection/ sur ili of 
tests, and design of underground defense structures. 


TI's talent for these “down-to-earth” studies dates back 
to 1930, when the company was formed as Geophysical 
Service Inc. — the first independent company to perform 
seismic reflection surveys for oil exploration. And from 


GSI’s early need for sensitive electronic seismic equipment 
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came an engineering skill that was applied to Anti-Sub- 
marine Warfare systems in World War II, and resulted in 
TI's Apparatus division becoming a leading ASW equip- 
ment manufacturer today. This rare blending of earth 
sciences knowledge and experience with skills in systems 
and components gives TI a unique capability in space-age 
technology. 


The founding company (GSI) is now the geophysical 
exploration arm of TI's Geosciences & Instrumentation 
division, and these capabilities are extended to govern- 
mental agencies by the Geosciences department. Designing 
and manufacturing advanced instruments and systems in 
wide use by these parties and others is the Instrumenta- 
tion Product group, with such products as WORDEN* 
gravity meters, EXPLORER* seismograph systems and 
seisMAC® seismic computers. 


For more about geosciences at TI, write to Central Mer- 
chandising for Bulletin DM-101. 


*Trademark of Texas Instruments 
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Now in the fourth year of operation, The Franklin Institute Com- 
puting Center has shared its technical know-how, modern electronic 
equipment and experienced personnel with hundreds of progressive 
industries and government agencies across the nation. 


We offer the services of creative people, highly skilled in their respec- 
tive fields and ably trained in the application of these skills to the ever 
expanding area of electronic computers and data processing systems. 


This staff is now available for analysis, system design, programming 
or coding of projects of unlimited scope or context. Input to our large 
scale computer and peripheral equipment is acceptable in any form. 
Results are provided on cards, plastic or metallic tape, and in com- 
pletely edited printed copy. Our extensive library of automatic coding, 
engineering, data processing and mathematical routines is available to 
all users, and machine time is provided with or without the services 
of programming personnel. 


THE FRANKLIN INSTITUTE 
Computing Center 


20th Street and Benjamin Franklin Parkway 
Philadelphia 3, Pa. LOcust 4-3600 
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